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SECTION  I 


INTRODUCTION 


The  results  of  the  NTS-2  independent  stability  and  control  analysis  will  be 
presented  in  the  following  order: 

o Summary  of  results  and  recommendations 

o Model  development 

o Simulation  and  supplementary  analysis 
o Conclusions 

Since  much  of  the  modeling  had  been  developed  under  the  NDS  study,  only 
those  model  elements  which  specifically  affect  NTS-2  results  will  be 
described. 


SECTION  II 


SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS 

NTS-2  verification  results  can  be  summarized  as  follows. 

Nutation  Damper  Performance 

Results:  Analysis  of  nutation  damper  performance  at  the  post  separation, 
post  AKM  burn,  and  post  AKM  jettison  mission  time  points  predicts  the 
following  ranges  of  equivalent  time  constants. 

• Post  separation,  300-700  seconds 

• Post  AKM  burn,  120-850  seconds 

• Post  AKM  jettison,  1000-5000  seconds 

The  ball  was  observed  to  hit  the  end  of  the  tube,  particularly  at  the  post 
separation  condition,  but  with  no  deleterious  effects  on  nutation  damping. 

Simulation  of  the  AKM  jettison  with  time  varying  inertias  and  worst  case 
disturbance  torques  indicated  an  insignificant  nutation  angle  disturbance. 

Recommendations:  None 

On  Orbit  Performance 

Results:  Large  stability  margins  were  predicted  from  linear  analysis  for 
the  effects  of  solar  torques,  magnetic  torques,  yaw  sun  sensor  errors  and 
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wheel  misalignments.  Simulation  of  satellite  performance  at  four  sun  line  to 
orbit  plane  angles  (0°,  12°,  45°,  86°)  under  nominal  and  worst  case  condi- 
tions uncovered  no  instabilities.  Wheel  failures,  if  compensated  for  within 
six  orbits,  will  cause  no  problems. 

Recommendations:  None 

Deployment  and  Acquisition 

Results:  If  reaction  wheels  are  turned  on  at  the  start  of  rod  deployment, 
earth/sun  acquisition  is  possible  without  earth/sun  sensor  saturation  for 
the  following  initial  body  rates 

p = . 007°/ sec  + h 

q = . 007°/sec  + h 

r = . 01°/sec  + y 

s 

where  k represents  earth  rate  and  y represents  sun  tracking  rate. 

s 

Recommendations:  Use  results  shown  above  as  guidelines  for  deployment 
procedures. 

Rod  and  Panel  Flexibility  Effects 

Results:  Large  gain  and  phase  margins  in  the  control  system  should  handle 
any  unpredicted  rod  flexibility  effects.  Simulation  of  rod  and  panel  flexibility 
dynamics  with  complete  satellite  model  uncovered  no  instabilities. 
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Recommendations:  None 


SECTION  III 


MODEL  DEVELOPMENT 

III -A.  REFERENCE  SYSTEM  DESCRIPTION 

Five  coordinate  systems  were  employed  in  the  analysis.  The  inertial  coor- 
dinate system  (Xj,  Yj,  Zj)  has  its  origin  at  the  center  of  the  earth.  points 
north,  Xj  is  positive  in  the  direction  of  vernal  equinox,  and  Yj  makes  up  a 
right  handed  triad.  The  local  vertical  coordinate  system  (X  Y Z ) 
has  its  origin  at  the  center  of  mass  of  the  space  vehicle.  Z^y  points  at  the 
earth's  center,  X^y  is  in  the  direction  of  orbital  velocity,  and  Y^v  is 
normal  to  the  orbit  plane.  The  inertial  and  local  vertical  coordinate  systems 
are  shown  in  Figure  1. 

The  body  coordinate  system  (Xg,  Yg,  Zg)  has  its  origin  at  the  center  of  mass 
of  the  space  vehicle.  Z points  in  the  direction  of  the  vehicle  antennae  or, 

D 

equivalently,  in  the  direction  of  the  gravity  gradient  rod  with  the  smaller 
tip  mass,  Y„  is 

D 

right  handed  triad.  This  axis  system  was  employed  with  the  NDS  analysis 
and  is  consistent  with  Rockwell's  definition.  A panel  coordinate  system 
(Xp,  Yp.Zp)  was  also  employed.  Z^  is  normal  to  the  panel  surface,  Y^  is 
along  the  solar  panel  support  boom,  and  X completes  the  right  handed 
triad.  Panel  angle  is  defined  to  be  zero  when  the  panel  coordinate  system 
is  aligned  with  the  body  coordinate  system.  The  body  coordinate  system  and 
panel  coordinate  system  are  pictured  in  Figure  2. 
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along  the  solar  panel  support  boom,  and  X_  completes  a 


Figure  1.  Inertial  and  Local  Vertical  Coordinate  System 
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Figure  2.  Body  and  Panel  Coordinate  System 


The  vehicle  data  provided  was  referenced  to  the  NRL  coordinate  system  which 
is  shown  in  Figure  3.  The  terms  in  parenthesis  represent  the  equivalent 
Honeywell  body  coordinate  system  translated  to  the  NRL  origin.  Vehicle 
mass  property  data,  as  provided  by  NRL,  is  given  in  Appendix  A. 

j 

III-B.  ONE  AND  THREE  BODY  MODELS 

The  NTS-2  verification  analysis  has  utilized  two  space  vehicle  models,  which 
have  historically  come  to  be  known  as  the  one  body  and  three  body  models. 

As  the  title  suggests,  the  one  body  model  is  less  sophisticated  than  the  three 
body  model.  In  the  NTS-2  verification,  the  one  body  model  was  used  in  all 
simulations  requiring  orbital  analysis.  It  assumes  the  following: 

• Space  vehicle  is  a rigid  body 

• Mass  properties  are  constant  (except  in  deployment) 

• No  dynamic  interaction  between  various  components  of 
space  vehicle  (except  reaction  wheels  on  space  vehicle) 

There  were  a number  of  other  approximations  made  to  alleviate  cost  prob- 
lems in  orbital  simulation  analyses;  however,  these  do  not  affect  the  basic 
structure  of  the  one  body  model.  These  will  be  discussed  in  Section  III-C. 

The  basic  structure  of  the  one  body  model  remained  the  same  for  the  NTS-2 
verification  as  for  the  NDS  analysis.  For  reference  purposes  it  is  included 
in  Appendix  B. 

Because  of  the  addition  of  the  gravity  gradient  rods,  the  three  body  model 
was  modified  from  its  NDS  counterpart.  Besides  providing  interaction 


Figure  3.  NRL/NTS-2  Coordinate  System 


between  all  components  of  the  space  vehicle,  the  three  body  model  also 
contains  flexibility  effects  of  the  gravity  gradient  rods  and  solar  panels, 
time  varying  inertia  properties,  and  time  varying  center  of  mass  position. 

There  are  some  space  vehicle  characteristics  that  can  only  be  analyzed 
with  the  three  body  model.  Nutation  damper  performance,  and  of  course, 
flexibility  effects,  cannot  be  analyzed  with  the  one  body  model.  The 
equation  development  for  the  three  body  model  is  given  in  Appendix  C. 

. 

III-C.  COMPONENT  MODELS 

Descriptions  of  the  components  of  the  NTS-2  model  are  now  presented. 
Assumptions  and  approximations  may  have  been  made  to  ease  simulation 
computational  burdens  and  will  be  so  stated  when  appropriate. 

III-C -1.  Nutation  Damper 

The  same  type  of  ball  in  tube  nutation  damper  used  in  NDS  is  used  in  NTS-2. 

The  physical  characteristics  are  slightly  different.  Figure  4 is  a Rockwell 
description  of  the  nutation  damper  for  NDS.  The  numbers  in  parenthesis 
refer  to  the  NTS-2  model.  For  NTS-2,  there  is  a specification  that  the 
nutation  damper  damping  coefficient  be  between  .01  and  .015.  (NRL  Nuta- 
tion Damper  Assembly  Design  Specification,  Specification  No.  NDL-SD-265498. 
Ford  Aerospace  Corporation,  25  April  1975.)  This  agrees  well  with  our 
predicted  estimate  of  performance  of  the  NDS  damper  which  has  a predicted 
damping  coefficient  of  .017  ± 30  percent.  A detailed  description  of  the 
nutation  damper  model  used  for  both  the  NDS  and  NTS-2  analysis  is  given  in 
Appendix  D.  All  pertinent  data  has  been  modified  to  reflect  NTS-2  charac- 
teristics. As  will  be  discussed  under  nutation  damper  performance  results, 
the  location  of  the  nutation  damper  does  have  a significant  effect  on  performance. 
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Figure  4.  NDS  and  NTS-2  Nutation  Damper  Characteristics 


III-C-2.  Attitude  Control  Electronics  System  Model 


Models  for  the  Attitude  Control  Electronics  System  (ACE)  were  constructed 
based  on  data  shown  in  Figures  5 and  6.  Two  versions  of  the  ACE  system 
were  constructed.  One  was  a complete  representation  of  the  ACE  system 
shown  in  Figures  5 and  6.  This  system  is  described  in  Appendix  E.  The 
system  was  used  when  the  effects  of  panel  and  rod  flexibility  were  studied. 

An  approximated  version  of  the  actual  system  was  also  constructed  to  pro- 
vide the  capability  for  making  simulation  runs  over  a number  of  orbits. 

This  was  done  to  avoid  the  need  for  the  small  integration  step-size  that 
would  be  required  if  the  actual  system  was  used.  The  rule  of  thumb  that 
we  use  in  choosing  step  sizes  for  the  integration  routine  is  to  select  the 
step-size  as  a function  of  the  highest  frequency  we  wish  to  simulate.  If 
there  are  poorly  damped  oscillatory  roots,  the  step -size  should  be  no  greater 

than  .1/u)  where  cd  is  the  largest  frequency  in  radians/sec.  The 
max  max 

largest  frequency  may  be  determined  by  looking  at  the  eigenvalues,  which 
in  the  closed  loop  case  may  be  represented,  as  follows: 

Source  Real  Freq  Damp 

Rigid  Body  Pitch  -.002 

-.000031 

Rigid  Body  Roll / Yaw  -.00187 

-.00187 

.000147  .1 

Each  Reaction  Wheel  .8  .62 

Each  Structural  Filter  -.033 
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Figure  5.  ACE  Control  System- -Spacetac  Print 


•e  6.  ACE  System — Spacetac  Print 


The  required  integration  steps  then  would  be  about  50  seconds  for  simulating 
rigid  body  dynamics,  3.0  seconds  if  structural  filters  were  included,  and 
. 125  seconds  if  the  actual  wheel  loops  were  included.  Hence  if  the  wheel 
loop  could  be  approximated  (and  the  structural  filters  ignored  because  they 
are  unnecessary  if  flexibility  is  not  modeled),  the  step  size  could  be  increased 
400  times  over  what  is  required  with  the  actual  control  system.  The  approx- 
imation was  derived  based  on  the  actual  intent  of  the  ACE  system.  Consider 
the  relationship,  as  shown  in  Figure  5 between  wheel  1 and  roll  attitude  error. 
This  may  be  represented  by  Figure  7. 

Since  the  wheel  loop  is  400  times  faster  than  the  dynamics  of  the  rigid  body 
and  since  we  neglect  the  structural  filters.  Figure  7 may  be  reduced  to 

W!  = K0  * 


COMPENSATION 

AND 

WHEEL 

J 

BLENDING 

DYNAMICS 

LOGIC 

(ROLL  ANGLE) 


(WHEEL  SPEED) 


Figure  7.  Roll  Angle  to  Wheel  1 Relationship 
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where  is  the  equivalent  transmission  gain.  However,  since  the  integra- 
tion is  with  respect  to  moment  equations,  we  utilized  the  relationship  between 
wheel  acceleration  and  wheel  torque  or 

Twl  = wi 


where  is  the  torque  wheel  1 produces,  and  1^  is  Wheel  1 moment  of 

inertia. 


Taking  the  derivative  of  equation- -and  substituting  in  equation — we  obtain 

Twl  = % Iwl  0 (2) 

where  1^.  is  the  loop  gain  for  that  wheel.  This  describes  the  basis  for 
the  approximated  version  of  the  ACE  system.  Block  diagrams  of  the  modified 
version  are  given  in  Figures  8 and  9.  As  shown,  there  were  several  other 
modifications  necessary  to  handle  the  non-linearities  in  the  ACE  system. 
These  will  be  discussed  as  follows: 


First,  the  effect  of  the  relays  at  the  output  of  the  coordinate  transformation 
matrix  shown  in  Figure  5 were  included.  This  was  done  by  monitoring  the 
output  of  that  relay  in  the  actual  control  system.  If  the  output  was  zero, 
which  indicates  a zero  commanded  wheel  speed,  it  was  assumed  that  the 
wheel  loop  could  be  unloaded  in  less  than  the  lime  of  an  integration  step.  In 
fact,  the  time  constant  of  that  loop  would  allow  a 95  percent  reduction  in 
wheel  speed  in  six  seconds  This  is  much  les-  than  the  50  second  integra- 
tion steps.  The  procedure  Emulated  when  tht  output  of  the  relay  was  zero 
was  to  command  a torque  which  would  drive  the  wheel  speed  to  zero  within 
that  integration  step,  or 
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Figure  8.  Modified  Wheel  Control  System 
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Figure  9.  Modified  Wheel  Dynamics 
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where  aT  is  the  integration  step. 


(3) 


Second,  the  effect  of  rolling  friction  was  modeled  by  defining  a maximum 
torque  limit  as  a function  of  wheel  speed,  as  shown  in  Figure  9. 

Third,  the  saturation  of  the  combined  earth  sensor  at  ± 4.0°  was  modeled 
• • _ 
by  setting  the  0 and  0 to  zero  when  0 and  0 was  greater  than  4.0  . 


Fourth,  the  yaw  sun  sensor  rate  signal  was  obtained  by  numerically  differ- 
entiating the  output  of  the  yaw  sun  sensor  model. 
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III -C  -3 . Magnetic  Torques  Model 

Spacecraft  magnetic  torques  result  from  the  interaction  between  the  mag- 
netic properties  of  the  spacecraft  and  the  ambient  magnetic  field.  The  main 
sources  are: 

• Permanent  magnetism  in  the  spacecraft 

• Spacecraft  generated  current  loops 

• Magnetism  induced  by  external  field 

• Currents  induced  by  external  fields 

• Rotating  assemblies 

The  earth's  magnetic  field  can  be  represented  by  an  analytical  model  that 
defines  magnitude  and  direction  as  a function  of  location.  Four  mathematical 
models  have  been  developed  for  representing  the  earth's  magnetic  field  and 
are  listed  below  in  order  of  increasing  accuracy. 

• Spin  axis,  centered  dipole 

• Tilted,  centered  dipole 

• Quadripole  model 

• Spherical  harmonic  expansion 

NASA  SP-8018  states  that  the  simple,  less  ac  t urate  models,  spin  axis  and 
tilted,  may  be  used 

• For  conceptual  or  preliminary  design  studies 

• When  the  magnitude  of  the  spacecraft  dipole  moment  M is 
an  estimated  value 
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• When  long-term  effects  are  of  primary  importance,  i.e. 
average  torque  or  torque  impulse  is  over  several  orbits, 
days,  or  longer 

• When  the  orbit  radius  (distance  from  the  geocenter)  is  greater 
than  three  Earth  radii 

In  addition,  data  given  in  SP-8018  indicate  a 5 percent  or  less  deviation  in 
magnitude  and  direction  between  a tilted,  centered  dipole  model,  and  the 
.sophisticated  spherical  harmonic  expansion  model  at  four  Earth  radii  (NTS-2 
orbit  radius). 

Based  on  these  two  considerations,  it  was  decided  to  use  a tilted,  centered 
dipole  model,  to  represent  the  earth's  magnetic  field.  The  details  of  the 
model  are  given  in  Appendix  F. 

III-C-4.  Solar  Torque  Model 

Mathematical  models  were  constructed  for  estimating  the  solar  torques 
produced  by  the  solar  panels,  the  gravity  gradient  rods,  and  the  tip  masses. 
The  basic  approach  was  to  break  up  the  component  of  interest  into  an  arrange- 
ment of  flat  surfaces.  The  solar  force  on  a flat  surface  may  then  be  com- 
puted by 

F = - FnN  + FtT  (4) 

where 

2 

Fn  = [(1  + a p)  cos  v + -p(i  - o)  1 1 ’A  cos  y 
F^  = (1  - a p)  cos  y sin  y PA 
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= normal  component  of  force 
N = unit  vector  normal  to  surface 
= tangential  component  of  force 
T = unit  vector  tangent  to  surface 

-7  2 

P = solar  pressure  constant  = 1.0  x 10  lb/ ft 

A = area  of  the  flat  surface 

p = fraction  of  incident  radiation  reflected 

a = fraction  of  reflected  radiation  that  is  specular 

(Note:  defined  another  way,  specular  reflection  = a o,  diffuse 
retlection  = p(  1 - a),  and  absorbitivity  a (1  - p)  ). 

Solar  torques  are  then  computed  by  summing  the  force-moment  arm  products 
with  respect  to  the  center  of  mass.  The  solar  torques  are  computed  as  a 
function  of  vehicle  attitude  and  panel  position  each  integration  step.  The 
nominal  value  used  for  the  radiation  characteristics  of  the  panels,  rods,  and 
tip  masses  are  given  in  Table  1. 

Table  1.  Radiation  Constants 


P 

CT 

Panels 

.2 

0 

Rods 

. 85 

1.0 

Tip  Masses 

0 

0 

For  more  details  of  the  solar  torque  model,  refer  to  Appendix  G. 
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Ill -C -5 . Gravity  Gradient  Hod  Thermal  Model 


The  effect  of  thermal  heating  on  the  gravity  gradient  rods  was  analyzed  using 
the  computer  program  MIT  AS  which  is  a software  package  provided  by  Control 
Data  Corporation.  MITAS  is  a general  three  dimensional  finite  difference 
heat  transfer  program  developed  to  mathematically  model  thermal  character- 
istics. The  approach  is  based  on  dividing  the  physical  model  into  an  R-C 
electrical  network  mesh  analog.  Either  explicit  or  implicit  finite  difference 
methods  can  be  used  for  steady  state  or  transient  solutions.  MITAS  can 
handle  all  three  types  of  heat  transfer,  conduction,  convection,  and  radiation. 

Deatils  of  the  model  can  be  found  in  Appendix  H.  The  major  results  may  be 
summarized  as  follows: 

• For  nominal  values  of  absorptivity  and  emissivity  (a  - 0.  15, 
e = 0.05)  the  predicted  radius  of  curvature  of  the  gravity 
gradient  rod,  assuming  the  sun  line  is  perpendicular  to  the 
rod,  is  2300  ft. 

• For  worst  case  values  of  absorptivity  and  emissivity 
(<*  = 0.30,  e = 0.1)  the  predicted  radius  of  curvature  is 
1100  ft. 

• Radius  of  curvature  varies  as  a esc  function  of  incidence 
angle. 

• Dynamics  of  thermal  deflection  have  time  constant  on  the 


III-C-6.  Gravity  Gradient  Rod  Flexibility  Model 
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A general  two-element  discrete  coordinate  structural  math  model  was  devel- 
oped for  simulation  of  a flexible  gravity  gradient  rod  for  NTS-2  satellite 
control  studies.  The  model  has  been  verified  with  STARDYNE  finite  element 
static  and  dynamic  check  cases.  It  also  allows  for  cross -axis  coupling 
effects  due  to  small  rod  curvatures  from  thermal  deflections  or  manufac- 
turing tolerances  on  deployed  rod  straightness. 

A series  of  multiple  beam  models  were  analyzed  for  dynamic  responses  to 
solar  array  drive  pulses.  These  studies  are  detailed  in  Appendix  I.  Based 
on  these  studies,  it  is  apparent  that  using  a single  flexible  beam  element  with 
only  one  node  at  the  tip  mass  gives  considerable  error  in  open-loop  transient 
response  simulations.  Using  four  elastic  beams  and  five  node  points  per  rod 
(one  attached  to  the  spacecraft),  realistic  results  are  obtained,  allowing  for 
realistic  rod  bending  modes.  Using  two  beams  and  three  nodes  per  rod  (one 
attached),  almost  identical  dynamic  results  are  obtained,  as  shown  in 
Appendix  I. 

Transient  responses  (open-loop)  are  presented  in  Figure  10  for  one  tip  mass, 
and  in  Figure  11  for  the  spacecraft  rotation  about  the  solar  array  drive  axis. 

Based  on  these  results,  a discrete  coordinate  math  model  for  the  gravity 
gradient  rods  was  developed.  Rod  deflections  at  the  tip  mass  and  at  a mid- 
point node  are  defined  in  terms  of  the  rod  stiffness  properties  (El  and  JG), 
rod  geometry  (2 £ = total  rod  length),  and  the  applied  loads.  For  deflection 
in  one  plane,  the  model  is  simply  a two  by  two  beam  influence  coefficient 
matrix.  For  cross -axis  coupling  due  to  rod  curvature,  the  initial  mass  point 
displacements  from  the  rod  undeformed  centerline  must  be  specified. 
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re  10.  NTS-2  Satellite  Tip  Mass  Acceleration  vs.  Time 


A<A 


Figure  11.  NTS-2  Satellite  Rotation  vs.  Time  About  Solar  Array  Drive  Axis 


Straight  Beam  Model --For  a straight  beam  model,  as  shown  in  Figure  12,  the 
elastic  deflection  of  nodes  A and  B are  as  follows: 


Figure  12.  Straight  Beam  Model 


where 

E = Young's  modulus  for  the  rod  material,  psi 

4 

I = rod  section  bending  stiffness,  in 


I = rod  half-length,  inches 


x,  y = elastic  displacements  in  x and  y directions  respec- 
tively, inches 

F , F = panel  point  loads  in  x and  y directions  respectively,  lbs. 

x y 


The  above  equation  assumes  the  same  El  for  bending  about  both  axes. 


The  panel  point  loads  are  defined,  in  any  one  axis,  as  follows: 

F = F - m x 

x x applied  relative 

where  F , . , could  be  the  x component  of  external  loads  from  solar 
x applied 

pressures,  or  the  x direction  component  of  a differential  gravity  field. 


, = the  acceleration  of  the  mass  point  as 

relative 

if  it  were  rigidly  attached  to  the  space- 
craft. 


Basic  Equations  of  Motion--In  terms  of  a simple,  one  degree  of  freedom 

harmonic  oscillator  subject -to-base  excitation  + applied  forces: 

• • • • . 

MX  + CX  + KX  = F(t)  - MX. 

base 


(6) 


Rewriting  this  equation 


KX  = -MX  - CX  + F(t)  - MX, 
• • 

X = e 


base 

MX  CX  , F(t)  M*b  ase 


K 


K + K 


K 


where  a = 1/K. 
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For  two  degrees  of  freedom,  the  basic  matrix  equation  is 


mA  ° " XA  + *CAACAB'  XA  + KAAKAB  XA  = RHg  (?) 
° mB  - XB  -CBACBB-  *3  -KBA  KBB  - XB 


It  can  be  shown  that  the  solution  to  the  above  equation  is 


= E 


mA  ° | XA 


o mT 


CAA  CAB 
CBA  CBB 


applied 

^applied 


m . o K 

A T base 


+ o mT 


base 


where  [E]  = influence  coefficient  matrix  = [K] 
and  [T"|  = rigid  body  transformation  matrix. 

Masses  are  defined  as  follows: 


m. . 
tip 

= tip  mass  = (7.5  lbs/gc>  for  NTS-2 

mbeam 

= rod  mass  = (.9  lbs/gc> 

mA 

= rn.  +1/4  m, 

tip  beam 

mB 

= 1/2  m 

beam 

lump 

= 1/4  m,  to  attach  point  on  spacecraft 

beam 

Curved  Beam  Model — Appendix  I contains  the  details  of  the  curved  beam 
deflection  calculations  for  unit  loads.  The  beam  Z axis  is  considered  to  be 
initially  aligned  along  the  vehicle  body-fixed  Z axis  but  then  is  displaced  to 
an  offset  unstressed  position  with  components 


X direction  offset  at  node  A 
X direction  offset  at  node  B 
y direction  offset  at  node  A 
y direction  offset  at  node  B 

The  model  is  shown  in  Figure  13. 

*0A  Z 


X 

Figure  13.  Curved  Beam  Model 


Note  that  the  offsets  are  considered  positive  if  in  the  negative  direction. 
(If  +X  is  toward  the  sun,  the  thermal  deflection  is  away  from  the  sun. ) 

The  full  influence  coefficient  matrix  is  presented  in  Table  2. 
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Table  2.  Influence  Coefficient  Matrix  for  Two-Beam 
Gravity  Gradient  Rod 
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INITIAL  UNSTRESSED  NODE  POSITION  AWAY  FROM  +Y  AXIS 
SUBSCRIPT  A REFERS  TO  TIP  MASS  NODE 
SUBSCRIPT  B REFERS  TO  NODE  AT  MIDDLE  OF  ROD 


For  linear  analysis  with  a curved  or  straight  beam,  the  influence  coefficients 
are  computed  once.  Numerical  values  are  presented  in  Appendix  I.  For  non- 
linear analysis,  the  influence  matrix  can  be  periodically  updated  if  the  cross- 
axis coupling  terms  are  updated  as  follows: 

Define  new  variables 

X ^ (current  value) 

X^g  = XqB  + Xg  (current  value) 

using  the  updated  variables  in  the  Table  2 equations. 

Damping --In  the  normal  method  of  calculating  structural  responses,  the 
damping  is  often  assumed  to  be  equivalent  viscous  damping,  for  lightly 
damped  systems.  The  modal  damping  is  conveniently  expressed  as  2 Q W , 
where 

Q = c/Cc,  the  damping  ratio  (usually  < . 10) 

W'  = the  mode's  frequency,  radians /sec 
Cc  = critical  damping  ratio 


If  the  modes,  frequencies,  and  generalized  weights  for  an  elastic  structure 
are  known,  it  can  be  shown  that  the  discrete  damping  matrix  can  be  calcu- 
lated using  the  estimated  modal  damping  values  (see  Appendix  I).  The  cal- 
culated damping  matrix  is 


(8) 


where 


C = discrete  damping  matrix 
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= square  matrix  of  mode  shapes 
= transpose  of  0 

= damping  ratio,  as  a fraction  of  critical 
viscous  damping 


S.  _ 


denotes  a diagonal  matrix 


= equivalent  or  generalized  mass  for  each 


mode  = 


= [ 0T  ] M 0 


In  accordance  with  previous  equations,  the  2 by  2 damping  matrix,  as  calcu- 
lated in  Appendix  I is 


C C 

'“AA  AB 

C C 

BA  BB 


201.06  -.75 


-.75  12.17 


x 10 


for  Q (or  g)  = . 005 


This  corresponds  to  1/2  of  1 percent  critical  damping,  which  is  an  estimate 
for  a gravity  gradient  rod.  The  full  damping  matrix  is  given  in  Table  3 in  a 
form  to  match  the  Table  2 influence  matrix. 


Table  3.  Damping  Matrix  for  Two  Beam  Gravity  Gradient  Rod 
(For  1/2  of  1 percent  Critical  Damping) 

fro  A 1 f 201.06  -.75  0 0 I f x 


-.75  12.17 


XB  (10> 


201.06  -.75 


-.75  12.17 
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where 

FD_  = damping  force  in  the  ith  direction,  at  the  jth  node,  lbs 

C = a conversion  factor  to  allow  for  possible  scaling  of 

above  damping  and/or  to  account  for  conversion  from 
inch  units  to  feet  units, 

x.,y.  = panel  point  velocities  at  the  jth  node  due  to  elastic 
^ ^ deflections,  inches /sec. 


1 I 

III-C-7.  Solar  Panel  Flexibility  Model 

A linear  structural  math  model  will  provide  reasonable  dynamic  simulations 
of  flexible  solar  array  panels  for  NTS-2  satellite  control  system  analysis. 

Bending  deflections  of  the  solar  arrays  are  obtained  from  the  discrete  mass, 
damping,  and  stiffness  matrices  described  in  Appendix  J.  This  approach  is 
similar  to  the  method  used  previously  to  model  the  flexible  gravity  gradient 
rods.  Separate  matrix  equations  are  given  in  Tables  4 and  5 for  each  array, 
each  consisting  of  three  panels.  Part  of  each  array's  mass  is  lumped  at 
the  center  of  the  outermost  panel.  The  balance  of  the  array  mass  is  added 
to  the  rigid  mass  of  the  spacecraft. 

The  torsional  flexure  model  cannot  be  simulated  very  well  with  a lumped 
mass  model  unless  a total  of  7 torsional  masses  are  used,  one  for  each 
panel  plus  one  for  the  drive  motor.  Therefore,  a modal  solution  is  de- 
veloped whereby  only  three  generalized  coordinate  responses  must  be  solved 
instead  of  seven  discrete  responses. 

The  preceeding  math  models  were  all  developed  on  the  basis  of  finite  element 
models  using  the  STARDYNE  structural  analysis  program. 
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Location  of  Drive  Motor--The  bending  equations  of  motion  are  defined  with 
respect  to  motion  of  the  solar  array  drive  motor,  which  is  considered 
centered  at  10.3  inches  from  the  CL  of  the  spacecraft,  located  on  the  +Y 
axis  (Honeywell  coord).  As  a result,  the  +Y  solar  array  has  short  drive 
tube  L =»  15,  inches,  while  the  -Y  solar  array  has  a drive  tube  about  35.5 
inches  long.  A sketch  provided  by  NRL  is  shown  in  Figure  14.  As  a result, 
the  bending  stiffnesses  for  the  two  solar  arrays  are  slightly  different. 

Bending  Stiffness,  Mass,  and  Damping  for  Solar  Arrays --Finite  element 
math  models  of  various  levels  of  complexity  for  the  solar  arrays  were  run 
on  STARDYNE,  as  reported  in  Appendix  J.  A model  with  a single  lumped 
mass  at  the  center  of  the  outer  panel  was  found  to  give  almost  the  same 
transient  responses  to  a step  yaw  torque  input  as  more  complex  models. 

The  general  solar  array  model  is  shown  in  Appendix  J.  Basically,  each 
solar  panel  has  a node  at  each  hinge  point  and  a node  at  its  center.  Addi- 
tional nodes  are  used  at  the  interfaces  with  a small  11  inch  long  honeycomb 
interface  panel  and  at  the  drive  tube  connections.  A total  of  32  nodes  was 
used  to  describe  the  flexible  solar  array  system. 

The  dynamically  equivalent  solar  array  bending  model  is  shown  in  Figure  15. 
The  corresponding  discrete  stiffness,  mass,  and  damping  matrices  are 
presented  in  Tables  4 and  5.  Note  that  the  x solar  panel  coordinates  are 
defined  to  be  normal  to  the  plane  of  the  solar  panels.  As  the  panels  are 
rotated  about  the  drive  axis,  the  local  x and  z panel  coordinates  also  rotate. 

The  slight  amount  of  stiffness  coupling  between  the  local  x and  z panel  deflec- 
tions are  attributed  to  a slight  asymmetry  in  the  hinges  between  the  solar 
panels. 
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Figure  14.  Solar  Array  Drive  Motor 


Figure  15.  Solar  Array  Bending  Flexure  Model 
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Table  4.  -Y  Solar  Array  (Bending  Only) 


STIFFNESS  MATRIX 
(units  = lbs /ft) 


For  motion  relative  to 
drive  motor 


[K]{x] « 


.0085617  -.004860  x 


•.004860  .630218  z 


where 

= motion  of  center  of  outer  panel, 
normal  to  plane  of  panel 

z^  = motion  of  center  of  outer  panel, 
in  plane  of  the  panel 


MASS  MATRIX 

(units  = slugs  (lb-sec2/ft)  ) 


[ M 3 f x 1 


.42176  0 " “x 

0 .52608  z 


DAMPING  MATRIX  q = .005  (1/2  of  1%  critical  damping) 


[C]  fx]  = 


.0006009 


.005758  z 


Units  of  [ C ] are 


lb-sec 
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Table  5.  +Y  Solar  Array  (Bending  Only) 


STIFFNESS  MATRIX 
(units  = lbs /ft) 


For  motion  relative  to 
drive  motor 


C K] {x  3 


.008670  -.010705 


-.010705  1.37791 


X2 

- L z9  . 


MASS  MATRIX  g 

(units  = slugs  (lb-sec  /ft)  ) 


C M ] (x  } - 


.42176 


.52608 


■ r xx  ' 
• • 

l zo 


DAMPING  MATRIX 
(units  = lb-sec /ft) 


[C]{ 


’ .0006047  0 If  x ‘ 

* i z 

x]  = 

0 . 0085 14  J z . 


The  damping  matrices  are  calculated  on  the  basis  of  the  uncoupled  mass 
and  stiffness  terms  for  each  direction.  For  a single  DOF  system,  with 


mass  m and  stiffness  k,  the  critical  damping  coefficient  Cc 

C = 2\fKm 
c 


is  defined  as 

(11) 


In  terms  of  the  equivalent  viscous  damping  ratio  £ , the  discrete  damping 
coefficient  is 


C = 2£Cc  (12) 

C = 2 Q\J  Km 


For  the  math  model  shown  in  Figure  16,  the  STARDYNE  model  was  subjected 
to  a step  yaw  torque  pulse  of  .04  ft -lbs  for  3 seconds,  after  which  time  the 
yaw  torque  goes  to  zero  instantly.  The  response  of  the  lumped  outer  panel 
mass  is  shown  in  Figure  17  in  inch  units.  Table  6 presents  summary  of  the 
peak  x direction  displacements  of  node  2 and  the  time  at  which  they  occur. 

It  demonstrates  that  the  equivalent  single  bending  D.  O.  F.  model  is  reason- 
able. Mode  shapes  for  bending  and  torsion  modes  are  presented  in  Appendix  J. 


Solar  Array  Torsional  Dynamic  Model- -A  simplified  torsional  dynamic  model 
for  the  solar  arrays  was  attempted,  using  a torsional  inertia,  lumped  at  each 
outer  and  inner  panel,  plus  a coordinate  at  the  drive  motor  location. 


The  first  symmetric  solar  array  torsional  mode  is  calculated  @ about  1.05 
Hz  (see  Appendix  J).  This  mode  has  the  inner  panels  plus  drive  motor  out 
of  phase  with  the  outer  panels.  This  mode  is  simulated  fairly  well  by  the 
simplified  model,  which  has  the  same  primary  torsional  mode  at  1.088  Hz. 
However,  the  second  symmetric  torsional  mode  at  1.89  Hz  is  missing  from 
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BEAMS  7,  8,  9: 

1.25  O.D. 

.049  WALL  THICK 
ALUMINUM  TUBE 


DRIVE  MOTOR 
ATTACH  TO 
S/CRAFT 


Figure  16.  STARDYNE  Math  Model  for 
Flexible  Solar  Arrays 


■ 
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Figure  17.  Gravity  Gradient  Rod  and 
Spacecraft  Definition 
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Table  6.  NTS-2  Solar  Arrays 


the  simplified  model  because  it  has  no  inertia  at  the  middle  panels  (nodes  5 
and  28).  The  middle  panels  have  the  maximum  amplitude  in  this  higher  mode. 

A series  of  DYNRE1  runs  were  made  in  the  complete  solar  array  model,  in 
which  the  effect  of  not  including  the  second  torsional  symmetric  mode  was 
compared  to  a run  with  all  significant  modes  considered.  Results  are  plotted 
in  Figures  18  through  23, 

Therefore  a modal  coordinate  approach  was  developed.  The  approach  is 
outlined  in  Appendix  J,  with  the  responses  of  the  generalized  coordinate 
computed  using  the  solution  to  the  basic  matrix  equation  for  2 degrees  of 
freedom.  If  the  discrete  torsional  response  of  the  drive  motor  or  the  center 
of  the  outer  panel  is  required,  perform  the  matrix  multiplication  indicated  in 
the  basic  matrix  equation  for  2 degrees  of  freedom. 

Solar  Array  Torsional  Dynamic  Model  Modifications--The  solar  array  tor- 
sional dynamic  model  was  modified  from  the  description  given  above  and  in 
Appendix  J.  The  modification  was  necessitated  by  a last  minute  change  in 
the  solar  array  hardware  which  significantly  increased  the  torsional  stiffness 
of  the  array.  A complete  description  of  the  modification  is  given  in  Appen- 
dix K.  The  biggest  impact  of  the  modification  is  in  the  frequencies  of  the 
first  two  symmetric  torsion  modes.  A comparison  between  the  original  and 
modified  mode  frequencies  is  given  in  Table  7. 

Also  shown  in  column  3 of  Table  7 is  the  effect  of  a grounded  drive  shaft  at 
the  drive  motor  versus  the  free  wheeling  model  frequencies  given  in  column 
2.  Although  the  shaft  is  not  completely  locked  or  grounded  at  the  drive 
motor,  there  are  torques  (Ball  brother  braking  torque,  rolling  friction) 
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Figure  18.  STAR  DYNE  Model  Node  2 
0Y  Response  Versus  Time 
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Figure  19.  Solar  Array  Torsion  Model  4 Pulse 
Input  Forcing  Function  (Normalized) 
(Amplitude  of  1.0  = .04  ft -lb) 
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Figure  20.  Solar  Array  Torsion  Model,  Node  2 = Outer  Panel 
Display,  Velocity,  Acceleration  and  Response  (in 
Radian  Units)  to  Four  Pulse  Excitation 
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Figure  23.  Response  of  -Z  Gravity  Rod  Tip  Mass  (Node  48) 

X Direction  Response  (Inch  Units)  Relative*  to 
Spacecraft  Attach  Point,  4 Pulses  on  Solar  Array 
Drive  Motor 

♦ NOTE:  No  rigid  body  contribution  for  relative  response. 
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Table  7.  Comparison  of  Original  and  Modified 
Torsional  Mode  Frequencies 


Mode 

Original  Model 
Frequency 
rad/sec  (Hz) 

Modified  Model 
Frequency 
rad/sec  (Hz) 

Grounded  Drive 
Motor  Model 
rad/sec  (Hz) 

1st  Symmetric 

6.88  (1.095) 

50.036  ( 7.96) 

16.28  (2.59) 

Torsion 

2nd  Symmetric 

11.925  (1.898) 

97.907  (15.53) 

56.77  (9.04) 

Torsion 

which  inhibit  motion  of  the  shaft  when  no  motor  torque  is  applied.  To  inves- 
tigate this  effect,  an  additional  solar  array  torsional  dynamic  model  was 
developed  which  included  the  drive  shaft  grounding  effects.  This  model 
differs  from  the  freewheeling  model  is  that  the  rigid  body  can  be  affected 
by  torsional  dynamics  of  the  array,  whereas  the  freewheeling  model  excites 
the  rigid  body  only  through  drive  motor  torques.  A complete  description 
of  the  additional  model  is  given  in  Appendix  L. 

III-C-8.  Gravity  Gradient  Rod  Deployment  Model 

During  rod  deployment  the  two  rods  are  extended  simultaneously  using  two 
dc  motors.  Each  rod,  however,  may  be  extended  or  retracted  independently. 
The  nominal  extension/ retraction  rate  is  1.25  in/sec.  Full  deployment  then 
would  require  approximately  9.6  minutes. 


Modeling--Deploying  the  gravity  gradient  rods  produces  a large  change  in 
the  moments  of  inertia  of  the  vehicle.  Hence,  the  gravity  gradient  rod 
deployment  model  consists  of  a time  varying  inertia  tensor  representation. 


The  inertia  tensor  for  the  total  system  may  be  represented  by 
* *RB  + ^panels  + *rodl  + *rod2 

The  first  two  terms,  the  rigid  body  inertia  tensor  and  solar  panel  inertia 
tensors  are  available  from  the  vehicle  mass  properties  data.  The  inertia 
tensor  of  the  rods  can  be  computed  based  on  the  theorem  of  parallel  axis. 


The  inertia  tensor  for  a hollow  right  circular  cylinder  about  its  principle 
axis  located  at  the  center  of  the  mass  of  the  cylinder  is  given  by 


*rods 


where 


mr  2 2 2 

Tf  tro  + ri  + 1 ’ 


m 


r r 2 2 2 

12  'ro  tri  +‘l 


mr.  2 2 

T tro  + ri  1 


(13) 


r = outside  radius  of  rod  = . 25  in 
o 

r.  = inside  radius  of  rod  = . 248  in 
mf  = mass  of  rod  = . 028  slugs 
j i = length  of  rod  = 60  ft  (fully  extended) 


The  inertia  tensor  for  a solid  shere  about  its  principal  axis  about  its  princi- 
pal axis  about  its  center  of  mass  is  given  by 
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*Tip 


mass' 


■=■  m r 
5 s s 


2 2 
=-  m r 
5 s s 


2 2 
■=■  m r 
5 s s 


(14) 


where 


m = mass  of  tip  mass  (m  = m = .233  slugs) 


r = radius  of  tip  mass  (r  =2  in;  r_  = 3.81  in) 

O * Q S _ 


The  principal  axis  of  the  tip  masses  can  be  chosen  so  as  to  be  parallel  to 
the  principal  axes  of  the  gravity  gradient  rods.  These  in  turn  can  be  con- 
sidered parallel  to  the  rigid  body  coordinate  system  which  are  approxi- 
mately vehicle  principal  axes. 


The  theorem  of  parallel  axis  can  now  be  used  to  calculate  the  contribution 
of  the  gravity  gradient  rods  to  the  vehicle  inertia  tensor  with  respect  to  a 
center  of  mass  coordinate  system.  The  total  system  inertia  is 


I 

xx 


I + I +1'  + m_  (r) 

xx  ....  . xx  , xx  , R. 

rigid  body  panels  rod^  1 

« 2 o 2 

+ I'  +m  (r  ) +1'  + (r  ) 

XXTip  mass i S1  XXrodj  2 

2 2 

+ 1'  + m (r  ) 

XXTip  mass2  S2 


(15) 
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where 


r = y component  plus  z component  of  distance  from  center  of 
mass  of  rigid  body  to  center  of  mass  of  rod  1 


r = y component  plus  z component  of  distance  from  center  of 
mass  of  rigid  body  to  center  of  mass  of  rod  2 


-1 


r = y component  plus  z component  of  distance  from  center  of 
mass  of  rigid  body  to  center  of  mass  of  tip  mass  1 


r = y component  plus  z component  of  distance  from  center  of 
mass  of  rigid  body  to  center  of  mass  of  tip  mass  2 


Similar  expressions  can  be  computed  for  I and  I 

yy  zz 


The  products  of  inertia  are  defined: 


I = I + I - m_,  (r  ^ r 

Xy  Xyrigid  body  xypanels  R1  x y 


1 ^ lv  < 2 2X  . ~ 2 ~ 2. 

“ mo  r ) - mr>  (r  r ) - rn  (r  r ) 


s,  x y 


R9  x y 


1 ~ J “2  “ J s2  X y 


I - I + I - m (r  r ) 

xz  xz  . , , xz  , R ' x z 

rigid  body  panels  1 


1 ~ 1.  ,2  2.  2 ~ 2. 
- m (r  r ) - m (r  r ) - m ( r r ) 
x z Rg  x z Sg  x z 


I = I + I - m„  (r  ^ r 

yz  yzrigid  body  yZpanels  R1  y z 


1 ~ 1.  ,2  2.  ,<v  2 ~ 2. 

- m (r  r ) - m (r  r ) - m (r  r ) 
Sj  y z R2  y z s2  y z 


where,  for  example. 


r^  , is  the  x component  of  the  vector  from  the  center  of  mass 
of  the  vehicle  to  the  center  of  mass  of  rod  1 


(16) 


J 
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The  center  of  mass  vectors  can  be  represented  as  follows: 


1 

r 


1 

P 


+ 


? 


1 


where 


1 

P 

1 

CT 


= vector  from  center  of  mass  of  rigid  body  to  NRL  origin 
= vector  from  NRL  origin  to  support  point  of  rod  1 


5 


1 


^3 


1 


vector  from  support  point  to  center  of  mass  of 
rod  1 


where 


l 

l 


1 

1 


= J ll  dt  Xj  £ = 60  ft 

max 

= constant  = rod  1 extension  rate  = 1.25  in/sec 


A similar  representation  can  be  given  for  rod  2. 


For  the  tip  masses 


~1  1 . 
r = r + 


0 

0 


+ l 
2 


The  masses  of  the  spheres  in  the  total  system  inertia  are  constant, 
ever,  the  masses  of  the  rods  change  as  they  are  extended. 


(17) 


How- 


54 


u 

0 


The  mass  of  the  rods  may  be  computed  as  follows: 


(18) 


m = YV 
1 

where 

Y = density  of  the  rod  material 

Y = volume  of  a hollow  cylinder 

Equation  m = yV  may  be  written 
rl 

mi  * 327TT  v'  I*-/  - ^ 

where  r^,  r.,  and  have  been  defined  previously  and 

Y = 690.  18  lbs/cu.  ft. 

Ill -C  -9 . AKM  Jettison  Disturbance  Torque  Model 

I 

Equipment  jettison  usually  produces  a small  transient  disturbance  torque 
that  has  a negligible  effect  on  the  long-term  performance  of  the  control 
system.  In  this  instance,  the  magnitude  of  the  disturbance  torque  effect  on 
spin  axis  attitude  is  of  prime  concern.  When  the  distance  between  the  jet- 
tisoned equipment  and  the  spacecraft's  mass  center  is  large,  the  line  of 
action  of  the  separation  force  should  pass  as  close  to  mass  center  as 
possible. 

In  the  present  case,  the  AKM  is  a spring-loaded  payload.  At  time  of  jet- 
tison, a retaining  band  is  released,  enabling  the  spring-loaded  payload  to 
be  ejected.  The  springs  do  not  expand  at  identical  rates  and  do  not  have  the 
exact  spring  constants. 
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In  addition,  it  cannot  be  assumed  that  all  springs  will  expand  exactly.  There- 
fore, there  will  be  a resulting  torque  acting  on  the  space  vehicle. 

Figure  24  illustrates  the  locations  of  the  four  springs  which  attach  the  AKM 
to  the  space  vehicle.  Important  data  is 


o 

Spring  constant: 

16  lbs/inch  ± tolerance 

o 

Uncompressed  length: 

1.5  inch  ± tolerance 

o 

Compressed  on  vehicle: 

. 75  inch  * tolerance 

o 

Solid  length: 

. 65  inch 

Figure  25  shows  the  positioning  of  the  springs  with  respect  to  the  X-Y 
coordinates  of  the  space  vehicle  assuming  no  AKM  misalignment.  Figure  26 
shows  the  same  positioning  assuming  the  AKM  misaligned  from  both  the  X 
and  Y axes.  Figure  27  illustrates  the  geometry  involved  in  calculating  the 
displacement  of  the  springs  from  the  space  vehicle  X-Y  axis,  assuming  no 
misalignment  (positioning  assuming  misalignment  can  be  easily  calculated 
from  these  base  calculations).  Table  8 contains  these  spring  force  locations. 

From 

F = kx 

the  force  created  by  a spring  expanding  .75  in.  with  k.  = 16  Ib/in  is  12  Lbs. 
The  resultant  force  of  the  four  springs  is  48  lbs.  The  weight  of  the  AKM 
is  64  lbs.  Its  mass  is  64/32.2  = 1.99  slugs.  Calculation  of  the  AKM 
acceleration  gives 

F = ma 
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Table  8.  Spring  Force  Displacement  and  Spring  Constant 


1 


48  = 1.99  a 

, 2 

a = 24. 12  ft/sec 

From 

d = *at2 

We  have 

.75/12  = %( 24. 1 2 )t 2 
t = . 072  sec 

and 

v = at  = 24.12  x .072  = 1.74  ft/sec 

The  resulting  velocity  of  the  AKM  is  calculated  to  be  approximately  1.74 
ft  /second. 

Torque  Computation- -Assuming  that  the  spring  constants  vary  by  no  more 
than  10  percent,  the  resulting  forces  can  be  calculated  and,  from  these 
forces,  the  disturbance  torques  can  be  easily  computed. 

Nine  cases  were  chosen  to  be  evaluated.  Case  1 involved  no  AKM  misalign- 
ment. Cases  2,  3,  4,  and  5 assumed  Y axis  misalignment  and  X axis  mis- 
alignment respectively.  Cases  6 through  9 assumed  both  X and  Y axis  mis- 
alignment. Table  8 lists  each  spring  constant  (including  a i 10  percent 
tolerance). 

Table  9 lists  the  computed  torques  about  the  X and  Y axis  of  the  space  vehicle. 
NOTE:  There  was  no  disturbance  about  the  Z axis. 
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Table  9.  Spring  Force  and  Resulting  Torque 


III-C-10.  Solar  Array  Drive  Model 


The  solar  array  drive  torque  model  is  based  on  data  generated  by  the  Ball 
Brothers  Research  Corporation.  The  data  is  in  the  form  of  handwritten 
notes  labeled  "Equations  for  Simulating  the  Rockwell  GPS  Solar  Array  Drive 
by  Glen  King  and  Lew  Martin."  It  is  dated  February  5,  1976.  The  pulse 
frequency  modulator  logic  is  as  described  in  a memo  to  Honeywell  from 
Samso  (P.  Talley,  October  22,  1975). 

General  Description — A functional  block  diagram  of  the  solar  array  posi- 
tioning loop  is  given  in  Figure  28.  Shown  is  a block  for  the  modification 
developed  by  GE  which  will  be  described  in  the  flexibility  analysis  section. 

Motor  Torques--As  given  in  the  Ball  Brothers  data  there  are  five  sources 
of  torque  associated  with  the  solar  array  drive  motor.  These  are: 

• Mechanical  torque 

• Coulomb  friction  torque 

• Cogging  torque 

• Magnetic  hysteresis  torque 

jjc 

• Flex  cable  torque 


M 


CL 


CG 


D 


FC 


In  addition.  Ball  Brothers  has  modeled  the  interaction  between  the  motor 
shaft  and  array  as  a torsional  spring  mass  damper  system  as  shown  in 
Figure  29. 


* 


This  term  relates  only  to  the  NTS-2  version  of  the  model. 


Figure  28.  Solar  Array  Drive  Functional  Block  Diagram 


Figure  29.  Torsional  Spring  Mass  Damper  System 


With  this  model,  motor  drive  shaft  angle  may  be  different  from  array  angle. 
The  dynamics  representing  this  interaction  are  given  by 


TS  *S  9S  " D(eA  " es’  ' K(6A  " es' 


«es-9A,+K(9s- VsIA«A 


where 


T = Total  torque  on  the  drive  shaft 
0g  = Shaft  angular  position 
0^  = Array  angular  position 
Ig  = Moment  of  inertia  of  shaft 
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1^  = Moment  of  inertia  of  array 

D = Equivalent  damping  coefficient 
K = Spring  constant 


The  total  torque  on  the  shaft  is 


rp  _ rip  rp  rTi  rri  rp 

s m cl  cg  d"afc 


(20) 


The  torque  producing  components  of  this  equation  have  been  modeled  by  Ball 
Brothers  from  experimental  data. 


They  are  (oz  - in) 

Motor  Torque 

tm*(25-‘-44  9s'<T7WTTt) 

= 8.48  (25  - 1.44  0 ) 

Coulomb  Friction  Torque 

T_,  = 30.0  For  > 0.0007  rad/sec 

LL  , b 

= 36.0  For  0C  < 0.0007  rad /sec  (during  pulse) 

= 0.  For  0g  < 0.0007  rad /sec  (no  pulse) 

Cogging  Torque 

Tcg  = [2.  3933  [sin  97  0g  + sin  (970g  + 0.803)] 

+ 0.3506  sin  (160g  + 1.57)]  {0.0385  + 0.964  sin2  (O.50g)] 

Magnetic  Hysteresis  Torque 

TDeJ*TDdt 


Td  = 14.6  [Td  SGN  0g  - 10.0]  0g 
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1 


Flex  Cable  Torque 


Tfc  = 4.0  sin  (0.9  0g)  (21) 

5 

Model  Simplification- -Values  for  the  spring  constant  (K  = 1.905X10  oz  - in/ 
rad)  and  equivalent  damping  coefficient  (D  = 352  oz-in-sec)  are  equivalent 
to  a torsional  mode  frequency  of  230.46  rad/sec  (36.68  Hz)  and  a damping 
coefficient  of  .213.  Preliminary  analysis  indicated  that  this  frequency  is 
much  higher  than  the  highest  significant  structural  mode-frequency  of  the 
panels  (less  than  20  rad/sec).  Thus  we  assumed  a rigid  connection  between 
the  motor  and  array.  A comparison  of  array  angular  response  for  the  motor 
drive  model  as  given  by  Ball  Brothers  and  the  rigid  approximation  is  given 
in  Figure  30. 

Array  Angle  Definition  and  Limits--NRL  defines  array  angle  to  be  zero  when 
the  array  is  in  the  yz  plane  (HI  coordinate  system)  and  with  the  solar  cell 
surface  pointing  in  the  +x  direction.  Positive  rotation  is  about  +y.  For 
GPS,  zero  array  angle  is  panels  in  the  xy  plane,  solar  cell  surface  pointing 
in  the  +z  direction.  Thus  for  panel  angle  0, 

0Hl"flNRL+9°  (22) 

The  limits  in  rotation  in  the  NRL  system  is  ± 100°.  In  our  system,  there- 
fore, 

- 10  s ©nj  £ +190°  (23) 

Since  the  motor  drive  is  not  situated  at  the  center  of  the  drive  shaft,  our 
structural  analysis  indicates  that  this  is  an  average  value  for  the  two  arrays. 


Figure  30.  Solar  Array  Response 


Pulse  Frequency  Modulator- -The  pitch  sun  sensor  signal  which  is  propor- 
tional for  ± 30°  (degrees)  (i.e.,  input  = output)  is  scaled  at  7.2  p amp/degree. 
This  signal  rescaled  drives  an  integrator.  A test  level  is  established  such 
that  the  integration  of  a 10°  error  for  .5  seconds  or  a 1°  error  for  5 seconds 
will  reach  the  same  test  value. 

Using  the  integrator  and  a timer,  variable  pulse  rates  between  . 2 Hz  and 
2 Hz  are  generated  with  the  following  logic: 

• The  integrator  and  timer  are  initialized  each  time  a pulse 
is  generated  or  the  timer  reaches  5 seconds. 

• If  the  error  signal  is  > 10°  (i.e. , the  test  level  would  be 
reached  in  < .5  seconds),  the  pulse  is  inhibited  until  the 
timer  reached  .5  seconds.  (This  limits  the  pulse  rate  to 
2 Hz.) 

• If  the  error  signal  is  < 1°,  the  test  level  will  not  be  reached 
before  the  timer  reaches  5 seconds  and,  hence,  no  pulse  is 
generated.  This  forms  a ± 1°  deadband  around  null. 

The  pulse  torque  polarity  is  determined  by  the  pitch  sun  sensor  polarity  and 
is  directed  to  rotate  the  drive  such  that  the  error  signal  is  reduced. 

Ill -C -11.  Reaction  Wheel  Model 

The  block  diagram  shown  in  the  Spacetac  blueprint  given  in  Figure  5 includes 
a description  of  the  wheel  loop.  The  effect  of  coulomb  friction  is  not  explicit. 
However,  it  has  been  included  in  the  feed  forward  block  containing  the  first 
order  lag  term  KtAt  S+l).  A more  realistic  approach  is  shown  in  Figure  31. 
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— may  be  approximated  with  data  shown  in  Figure  32  which  was  taken  from 
Bendix  test  data*  or 


T T 

&T  _ 1 - 2 _ 2.1  in  oz 

bw  w - w2  -1330  rpm 


(25) 


which  assumes  linear  operation  over  the  total  operating  speed  of  the  wheel. 
The  time  constant  then  may  be  calculated  as 

J_  = -OH  ^ lbs  1 (2f 

tw  139. 28  rad/sec  .0031  slug  ft" 

= .0255  sec  * 


or 

t = 39.25  sec 
w 


The  lines  labeled  +20%  linearity  is  equivalent  to  a 33  second  time  constant 
and  the  -20%  linearity  is  equivalent  to  a 52  second  time  constant.  Thus  T 

w 

could  vary  between  30  and  50  seconds. 


III-C-12.  Combined  Earth  Sensor  Model 

A mathematical  model  of  the  combined  earth  sensor  which  includes  both  a 
horizon  crossing  indicator  for  use  in  the  spin  phase  and  %n  earth  radiation 
detector  for  on  orbit  generation  was  developed  for  the  NDS  verification.  The 
model  includes  a filter  with  a 1.4  second  time  constant  representing  thermo- 
pile image  processor  fade  lag.  The  sampling  strategy  for  computing  pitch 


*NRL  reaction  wheel  type  No.  3890003-1,  Bendix  Navigation  and  Control 
Division,  28  August  1974. 
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POWER  TORQUE 

(WATTS)  I (OZ-IN) 


Figure  32.  Bendix  Test  Data 


and  roll  errors  is  also  included.  A complete  description  is  given  in 
Appendix  M. 

For  the  NTS-2  analysis,  only  the  earth  radiation  detector  model  was  used. 
For  use  in  simulating  satellite  performance  over  many  orbits,  the  image 
fade  lag  and  the  sampling  strategy  were  removed  from  the  model.  The 
pitch  and  roll  signal  comparison  tests  were  retained  however. 

The  complete  earth  radiation  detector  model  was  used  in  the  simulation 
runs  which  investigated  rod  and  solar  array  flexibility  effects. 

in-C-13.  Yaw  Sun  Sensor  Model 

The  yaw  sun  sensor  is  mounted  on  the  middle  solar  panel.  Its  input -output 

characteristics  are  shown  in  Figure  33.  The  input  parameter,  rotation 

angle  from  boresight  (y  ),  is  calculated  from 

s 

, -i  V 

Y = -tan  yp 
S ~ 

zp 

where 

V = y component  of  the  sun  line  vector  in  panel  coordinates 
= z component  of  the  sun  line  vector  in  panel  coordinates 

The  output  has  been  transformed  from  p amps  to  degrees.  The  operating 
accuracy  of  the  sensor  at  25°C  is  bounded  by  the  dashed  lines.  There  is  an 
additional  scale  factor'1'  that  is  applied  to  the  output  to  simulate  temperature 

3$C 

Scale  factor  range  suggested  by  Rockwell. 
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Figure  33.  Yaw  Sun  Sensor  Characteristics 


variations.  Scale  factor  range  due  to  temperature  variations  were  assumed 
to  be  ±50%  from  the  nominal  value  of  unity. 

III-D.  MODEL  APPROXIMATIONS,  UNCERTAINTIES,  AND 
DESIGN  TOLERANCES 

Most  of  the  modeling  approximations,  level  of  uncertainty,  and  component 
design  or  alignment  tolerances  were  discussed  in  each  of  the  component 
model  development  sections.  This  section  will  present  a summary  of  these 
model  variation  items. 

The  major  approximations  were  necessitated  by  the  desire  to  simulate 
satellite  performance  over  a number  of  orbits.  Hence,  elements  of  the 
model  which  were  associated  with  relatively  fast  dynamics  were  either 
removed  or  modified.  The  approximations  made  for  these  orbital  nans 
consisted  of  the  following: 

• One  body  model  (i.e.  constant  inertia  except  for  deployment) 

• Modified  ACE  control  system 

• No  time  lag  or  sampling  strategy  in  earth  sensor  (most  runs 
made  with  actual  pitch  and  roll  feedback  with  4°  limit) 

• No  solar  array  drive  (perfect  panel  tracking) 

• No  panel  or  rod  flexibility 

The  three  body  model  which  was  used  to  investigate  flexibility  effects  con- 
tains no  approximations  to  the  models  as  originally  conceived.  We  feel 
that  it  is  as  complete  and  accurate  a model  as  possible.  Operating  ranges 
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of  pertinent  component  models  will  now  be  summarized.  Sources  of  design 
tolerances  or  alignment  specs  will  be  noted  when  applicable. 

Nutation  Damper 

• Damper  coefficient  0.01  - 0.015  (NRL  Spec) 

• Radius  of  curvature  13±1  ft  (NRL  Spec) 

Magnetic  Torques 

• Spacecraft  dipole  moment  500  - 2000  pole-cm 

Rod  Thermal  Curvature 

• Min  radius  of  curvature  1100  - 2300  ft 

Earth  Sensor  Errors 

• Gross  error  0.25°  (NRL) 

• Fine  errors  as  discussed  in  Appendix  M 

\ 

\ 

Yaw  Sun  Sensor  Errors 

• Gain  variation  ±50%  (Rockwell) 

• Envelope  as  discussed  previously 

Vehicle  Mass  Properties 

• CM  misalignment  Ax  = Ay  = . 05  in,  iz  * . 1 in 

2 

• Inertia  variation  I , I , I 2.2  slug  ft  (NRL) 

xx  yy  zz  b 

I 0.9  slug  ft^ 
xy  6 


1 
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Wheel  Misalignment 

• 1°  in  cant  and  yaw  (NRL) 

Panel  Center  of  Pressure  Shift 

• 0.5  inches 


Rod  Deployment  Rates 

• * 10  percent 

AKM  Jettison  Springs 

• - Spring  constant  ±10  percent 


11 
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• Nutation  damper  performance 

• On  orbit  performance 

• Deployment  and  acquisition 

• Rod  and  panel  flexibility  effects 

The  results  of  each  of  these  topics  will  now  be  analyzed. 

IV -A - 1 . Nutation  Damper  Performance 

Nutation  damper  performance  is  a function  primarily  of  satellite  inertia 
ratio,  nutation  damper  dissipation  coefficient,  and  location  of  the  nutation 
damper  with  respect  to  the  spin  axis.  As  shown  in  Appendix  D,  the  motion 
of  the  ball  in  the  tube  can  be  modeled  with  spring  damper  equations  where 
the  frequency  of  the  response  is  given  by 

%d  = 7 R % (28) 

where 

X = distance  of  nutation  damper  from  center  of  mass 

R = radius  of  curvature  of  the  tube 

a)  = spin  rate 
z 


Optimum  performance  of  the  nutation  damper  will  occur  where  the  nutation 
damper  frequency  is  tuned  to  the  nutation  frequency  which  is  given  by: 


Equating  ( ) and  ( ) 


5 l 
7 R 


tu 


'XX  ’ 


ou 


we  see  that  tuning  the  damper  is  independent  of  the  spin  rate.  Although 
nutation  damper  satellite  dynamics  are,  at  a minimum,  fifth  order,  nuta- 
tion performance  has  been  expressed  as  an  equivalent  first  order  time 
constant.  Figure  34  presents  plots  of  equivalent  time  constant  versus  1 
for  different  inertia  ratios.  The  time  constants  were  computed  uiing  the 
"energy  sink"  approximation  described  in  Appendix  D.  Also  shown  is  the 
location  of  the  NTS-2  damper  at  1.5  feet  rom  the  spin  axis.  At  that  loca- 
tion, the  nutation  damper  will  perform  optimally  for  an  inertia  ratio  of 


(DISTANCE  FROM  SPIN  AXIS,  FT) 

Figure  34.  Damper  Position 
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1.3.  For  inertia  ratios  greater  than  1.3,  performance  will  deteriorate 
much  more  rapidly  than  for  those  less  than  1.3. 

For  NTS-2,  nutation  damper  performance  was  investigated  at  three  mission 
time  points:  1)  post  separation,  2)  post  AKM  burn,  and  3)  post  AKM  jettison. 
This  represents  an  inertia  ratio  spread  of  1.17  at  post  separation  to  1.54 
after  the  AKM  has  been  jettisoned.  Nutation  damper  performance  was  sim- 
ulated at  these  three  time  points  under  nominal  and  worst  case  conditions 
for  spin  rates  of  90  RPM  (nominal),  100  RPM  (maximum  expected),  and  40 
RPM  (minimum  allowable).  The  worst  case  conditions  chosen  were: 

o Damper  dissipate  coefficient  -0.01 

2 

o Inertia  ratios  AI  =Al  = Al  =±2.2  slug  ft 

xx  yy  zz 

o Cm  misalignment  Ax  = Ay  = 0.1"  Az  = 0.2" 

The  sign  on  the  inertia  ratio  variation  was  chosen  such  that  the  time  con- 
stant variation  was  largest  as  estimated  from  Figure  34.  The  sign  on 
was  chosen  to  increase  the  distance  of  the  center  of  mass  of  the  nutation 
damper  from  the  center  of  mass  of  the  satellite. 


Nutation  damper  performance  in  terms  of  nutation  angle  time  history  for 
the  post  separation  configuration,  and  nominal  conditions  at  90  RPM  is 
shown  in  Figure  35.  As  shown,  the  estimated  equivalent  time  constant  is 
approximately  300  seconds.  (Note:  as  discussed  in  Appendix  D,  the  nuta- 
tion damper  model  assumed  perfect  elastic  collisions  when  the  ball  hit  the 
end  of  the  tube.  The  response  shown  in  Figure  35  was  for  a model  which 
assumed  the  other  extreme,  an  inelastic  collision.  A comparison  of  per- 
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Nutation  Damper  Run  Launch  90  RPM, 
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formance  will  be  discussed  shortly. ) The  results  of  nutation  damper  perfor- 
mance evaluation  at  the  remaining  mission  time  points  is  given  in  Figures  36 
through  38  and  in  Tables  10  through  12.  These  conditions  were  simulated 
at  the  Air  Force  Academy,  where  it  was  not  possible  to  generate  plots  of 
the  form  shown  in  Figure  35.  Hence  only  the  numerical  results  are  given. 

Equivalent  time  constants  estimated  using  two  other  techniques  are  also 

shown,  t is  the  predicted  time  constant  using  "energy  sink"  approxima- 
es 

tions.  t is  the  time  constant  computed  from  eigenvalues  generated  at 

At 

selected  points  in  the  simulation  via  a numerical  linearization/ root  fender 
procedure.  The  latter  values  were  used  as  a check  on  simulation  efficacy. 
Time  constants  calculated  using  the  energy  sink  approximations  generally 
predict  pessimistic  performance  particularly  at  the  larger  inertia  ratios. 

Since  there  is  no  defined  performance  specification  applied  to  the  nutation 
damper,  the  results  shown  in  Figures  36  through  38  and  Tables  10  through 
12  do  not  in  themselves  indicate  satisfactory  or  unsatisfactory  performance. 
That  evaluation  is  a function  of  mission  time  line  constraints.  To  aid  in 
making  those  evaluations,  a summary  of  the  results  is  as  follows: 

Time  Constant  (sec) 


Nom 

Worst  Case 

Post  Separation 

300 

700 

Post  AKM  Burn 

150 

850 

Post  AKM  Jettison 

1100 

5000 

In  addition,  the  ball  does  hit  the  end  of  the  tube  for  initial  nutation  angles  of 
3°  or  more  (post  separation).  This  does  not  appear  to  have  a significant 
effect  on  damping  whether  the  tube  collision  is  assumed  to  be  elastic  or 
inelastic.  The  elastic  collision  model,  which  was  used  in  all  simulation 
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TIME  CONSTANT  (SECONDS) 


TIME  CONSTANT  (SECONDS) 


Figure  37.  Damper  Performance,  Post  AKM  Burn 
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TIME  CONSTANT  (SECONDS) 


600  - 835  Yes 


Table  11.  Nutation  Damper  Performance — Post  AKM  Burn 


620  520  1062 


3000  - 8015 


runs,  was  verified  to  be  worst  case.  Although  the  difference  in  equivalent 
time  constant  between  the  elastic  and  inelastic  models  for  identical  runs 
was  small,  the  ball  hits  the  end  of  the  tube  twenty  seconds  longer  into  the 
run  with  the  elastic  model  than  with  the  inelastic  model. 

IV-A-2.  A KM  Jettison 

In  order  to  determine  what  effect  jettisoning  the  AKM  would  have  on  nutation 
angle,  the  jettison  was  simulated  using  the  worst  case  disturbance  torques 
selected  from  Table  9 . These  were 

ATx  = -1.834 

LTy  = -3. 178 

These  torques  were  applied  over  .075  seconds,  that  is,  the  time  it  takes 
the  springs  to  expand  to  their  uncompressed  length.  Time  varying  off 
nominal  inertia  properties  were  also  included  in  the  simulation.  The 
results  are  shown  in  Figure  39.  As  can  be  seen,  peak  nutation  angle  due 
to  the  AKM  jettison  is  approximately  .006°,  an  insignificant  amount.  It  is 
concluded  that  the  AKM  jettison  should  produce  no  undesirable  effects. 

IV-B-1.  On  Orbit  Stability  Analysis 

The  closed  loop  control  system  is  a combination  of  active  and  passive  con- 
trol. The  passive  control  moment  producers  are  the  pitch  and  roll  gravity 
gradient  torques.  The  active  control  moment  producers  are  four  reaction 
wheels  aligned  to  give  decoupled  three  axis  control.  As  discussed  in  the 
modeling  section,  the  ACE  system  commands  wheel  momentum  proportional 
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Figure  39 . Nutation  Damper  Run  AKM  Jettison 


to  attitude  error,  or  equivalently,  wheel  torque  proportional  to  attitude  rate. 
Commanded  wheel  momentum  may  be  represented  in  satellite  axes  as  follows: 


Hw 

a* 

ii 

0 

X 

= K 

6 

w 

y 

y 

Hw 

= K 

Y 

z 

and  similarly 


T ...  - K $ 

W x 

x 

TW  = Ky  6 

y 

Tw  = Kz* 

z 


where 


) = commanded  wheel  momentum  (torque)  in  roll  axis 

x x 

H^.  (T-^.  ) = commanded  wheel  momentum  (torque)  in  pitch  axis 

y y 

HW  <TW  ^ = commanded  wheel  momentum  (torque)  in  yaw  axis 
z z 

Kx  = roll  axis  loop  gain 

= pitch  axis  loop  gain 
K = yaw  axis  loop  gain 


In  inertial  coordinates,  total  torque  includes  the  effects  of  satellite 


orbital  rate  and  can  be  approximated  by 
= + 
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Substituting  equations  H^.  and  into  Tj,  we  get 


TI  = Kx  0 + TlKY 


ti  - K e 
y J 

Ti  3 Kz  ^ " ^Kx  * 


For  closed  loop  dynamics  this  equation  is  combined  with  the  solid  shere 
inertial  tensor  equation  to  form 

I 0 - K 0 - 4T]2  (I  - I ) 0 + (I  + I - I ) TI  Y - T|K  Y = 0 (33) 

xx  x ' zz  yy  zz  xx  yy  z 

I 0 - K 0 + 3t)2  (I  - I ) 0 = 0 
yy  y 1 xx  zz 

I Y-Ky  + T12(I  -I  ) Y - (I  -I  -I  )T10+11K  0 = 0 

zz  z 1 yy  xx  yy  xx  zz  1 x 

The  nominal  values  of  the  constant  in  the  above  equation  are  given  in  Table  13 
Table  13.  Satellite  Equation  Constants 


I 

1998  slug  ft2 

xy 

2 

i 

1914  slug  ft 

yy 

2 

i 

zz 

205  slug  ft 

K 

3. 75  ft-lbs/rad 

X 

K 

3.  75  ft-lbs/rad 

y 

K 

.4  ft-lbs/rad 

z 

TI 

.000145  rad /sec 
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The  nominal  eigenvalues  of  the  satellite  have  been  listed  previously.  Eigen- 
value migration  for  the  three  axes  is  shown  in  Figure  40  as  a function  of  loop 
gain.  Since  roll  and  yaw  are  coupled  the  loci  are  for  each  axis  with  the  other 
axis  closed  at  the  nominal  loop  gain. 


In  the  pitch  axis,  the  roots  migrate  from  the  open-wheel  loop,  neutrally- 
stable,  gravity-gradient,  oscillatory  pair  to  a pair  of  real  roots  at  -.002 
and  -.000031.  These  roots  result  from  a loop  gain  that  is  four  times  larger 
than  that  required  for  critical  damping. 


The  roll  and  yaw  loops  are  characterized  by  a low  frequency  poorly  damped 
oscillatory  pair  which  is  due  to  the  coupling  in  the  control  loop  through  the 
orbital  rate  term  as  shown  in  the  above  equation.  This  coupling  is  relied 
upon  for  stability.  This  can  be  seen  more  clearly  if  we  look  at  the  s plane 
representation  of  the  solid  shere  inertial  tensor  equation  in  matrix  form  as 
shown  in  Figure  41.  The  characteristic  equation  of  the  system  shown  in 
Figure  41  may  be  found  by  forming  the  determinant  of  G(s).  The  charac- 
teristic equation  will  have  the  form 


4 3 2 

s + a3s  + a2s  + alS  + aQ  = 0 


(34) 


The  importance  of  the  roll  yaw  control  coupling  can  be  seen  by  looking  at  af 


which  is 


^ 2 (Iyy  “ Izz)  \ 

[ '2  (Iyy  " Ixx)  ^ 

I 4. 

Ik  i 1 \ 
1 x 

i 411  ) 

r *«  j 

1 T 

i I 

1 

V zz 

(35) 


The  first  term  of  this  equation  is  independent  of  reaction  wheel  control. 
The  second  term  is  a direct  result  of  the  coupling  in  roll  and  yaw  wheel 
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PITCH  AXIS 


ROLL  AXIS 
(YAW  CLOSED) 

Kx  = 3.75 


YAW  AXIS 
(ROLL  CLOSED) 

<2  = 0.4 


Figure  40.  Root  Migration  as  a Function  of  Loop  Gain 


Figure  41.  Roll/Yaw  s Plane  Equations 


control  through  orbital  rate  as  given  before.  Since  I > I > I , if  K or 
b xx  yy  zz  x 

is  equal  to  zero,  then  a^  < 0 and  the  system  is  unstable. 

In  this  section,  we  have  attempted  to  briefly  describe  the  nature,  and  partic- 
ularly the  stability,  of  the  NTS-2  control  system.  The  reaction  wheel  loop 
gains  have  been  chosen  from  steady  state  attitude  error  considerations  under 
the  constraint  of  adequate  stability  margins  for  flexure  effect  and  uncertain- 
ties. The  final  gains  selected,  when  viewed  from  a root  locus  analysis, 
produce  low  frequency  roots  which,  at  first  glance,  were  uncomfortably 
close  to  the  imaginary  axis.  In  the  results  that  follow,  we  show  we  could  not 
destabilize  those  roots  under  reasonable  off  nominal  perturbation  conditions. 

Analyses  of  Destabilizing  Effects --As  stated  above,  the  poorly  damped  nature 
of  the  pair  of  oscillatory  roots  in  the  roll/ yaw  equations  caused  concern.  We 
tried  to  establish  whether  or  not  there  were  destabilizing  effects  that  could 
force  those  roots  into  the  right  half  plane.  The  destabilizing  effects  that  were 
investigated  were 


96 


•TpariiTiiirn 


o Yaw  sun  sensor  errors 

o Yaw  solar  torques 

o Magnetic  torques 

o Wheel  misalignments 

We  will  now  discuss  the  results  of  that  investigation. 

Yaw  Sun  Sensor  Errors --The  control  equations  described  previously 
assumed  that  the  yaw  sun  sensor  sensed  pure  yaw  rotations.  There  are  only 
two  points  in  each  orbit  where  this  is  the  case,  assuming  perfect  sun  track- 
ing on  array  pitch  control.  Since  the  yaw  sun  sensor  is  mounted  on  the  solar 
array  and  measures  yaw  as  the  angle  between  the  sun  line  and  the  array  rota- 
tion axis  in  the  plane  defined  by  the  array  rotation  axis  and  the  array  normal, 
the  sensed  yaw  angle  is  a function  of  array  angle  with  respect  to  the  satellite 
body.  The  sensor  senses  pure  yaw  only  when  the  array  normal  is  perpen- 
dicular to  the  satellite  yaw  axis.  A derivation  of  what  the  sensor  sees  for 
the  case  of  the  array  normal  being  tilted  45°  to  the  yaw  axis  now  follows. 

The  sensed  yaw  angle  is  given  by 


y = tan 
s 


. V 

-i  yp 


where 


V = y (array  axis  of  rotation)  component  of  sun  line  in  panel 
^ coordinates 


= E_  R + E R + E R 
21  x 22  y 23  z 
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Vyp  = z component  (array  normal)  of  sun  line  in  panel  coordinates 

= sin  0 V „ + cos  0 V _ 
p xB  p zB 


where 


E„  (i  = l,  3 j = l,  3)  = components  of  inertial  to  body  transformation 


R , R , R 
x y z 


= sun  line  in  inertial  coordinates 


= x component  of  sun  line  in  body  coordinates 
= z component  of  sun  line  in  body  coordinates 


To  force  0 = 45  R = -.  707,  R = 0,  R = -.707.  The  sun  line  in  panel 
p x y 2 

coordinates  then  becomes 


V =-(.707)(E21+E23) 

P 

V = . 707(  - . 707)  (E  + E - E„  - E„_) 

z 11  13  31  33 

P 


To  first  order. 


Y = tan 
s 


-i  % 


( + cos  y sin  0 sin  0 - sin  Y cos  0 + cos  0 sin  0 ) 

^ s . 707  (cos  0 cos  Y - sin  0 + cos  Y sin  0 cos  0 + sin  Y sin  0 

+ COS  0 COS  0) 

Assuming  small  angle  approximation 

' TTu  - *>  (40) 

It  is  seen  from  this  equation  that  when  the  solar  array  is  tilted  45°  with 
respect  to  the  satellite  yaw  axis,  the  sun  sensor  senses  an  equal  combination 
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of  roll  and  yaw.  At  other  array  angles  the  sensor  will  sense  different  pro- 
portions of  yaw  and  roll.  This  effect  was  investigated  by  modeling  the 
sensor  output  yg  as 

y = ay  + b0  (4 

s 

where  the  coefficients  a and  b represent  the  effect  of  array  rotation  angle. 
Substituting  this  equation  for  y in  the  solid  shere  inertial  tensor  equation 
and  ignoring  the  pitch  component,  we  obtain 

T = K 0 + ”T]K  (ay  + b0)  (4 

I X z 

X 

T = K (ay  + b0 ) - T]K  0 

L Z Z 

z 

The  G(s)  matrix  which  includes  this  effect  is  shown  in  Figure  42. 
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Figure  42.  G(s)  with  Roll/Yaw  Sensor  Feedback 


We  again  form  the  characteristic  equation  to  look  at  the  effect  of  the  yaw 
sensor  feedback  on  stability.  It  is  particularly  interesting  to  look  at  the 
a^  coefficient. 
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coupled  yaw  sensor  output  is  a factor  of  a gain  reduction. 

Although  the  same  cancellation  effect  does  not  occur  in  the  other  coefficients, 
it  is  the  rate  coefficient  that  has  the  most  impact  on  the  location  of  the  oscil- 
latory pair.  Table  14  shows  the  effect  of  varying  b on  the  location  of  the 
oscillatory  pair.  No  effect  can  be  seen. 

Table  14.  Effect  of  Varying  Yaw  Sensor  Output  on  Root  Migration 


a 

b 

Oscillatory 
Damping  Ratio 

1.0 

0 

. 1033 

.0001465 

1.0 

0.2 

. 1033 

.0001465 

1.0 

1.0 

. 1033 

.0001465 

L 


100 


0 


11 

u 

0 


0 

IJ 

li 

D 

c 

0 


( 
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The  effect  of  gain  reduction  can  be  seen  in  the  roct  locus  of  Figure  39.  The 
major  effects  is  on  the  location  of  the  real  root.  There  is  little  impact  on 
the  location  of  the  oscillatory  pair. 


Yaw  Solar  Torques --The  most  significant  effect  of  solar  pressure  is  a 
yaw  moment  which  is  produced  because  of  the  offset  of  gravity  gradient  rods 
from  the  center  of  mass  of  the  satellite.  This  is  shown  schematically  in 
Figure  43.  The  rods  are  normal  to  the  paper.  As  can  be  seen,  if  the  x axis 
of  the  satellite  is  perfectly  aligned  with  the  sun  line,  there  is  no  torque. 
However,  as  the  vehicle  yaws  there  is  a destabilizing  yaw  torque  created 
due  to  the  offset  of  the  rods  on  the  x axis.  This  can  be  modeled  for  stability 
analysis  purposes  as  follows: 


The  force  on  the  rods  is  approximately  given  by 
F = (1  + p)  PA 


(44) 


Figure  43.  Yaw  Solar  Torques 
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where 


p = solar  reflectivity  constant  = . 85 

“7  2 

p = solar  radiation  pressure  constant  = 1 x 10  lbs  /ft 

5 2 

A = area  of  rods  = 2 x 60  x = 5 ft 
The  torque  generated,  T , is  then 

Si 

Tst  = F X / sin  Y 

where  l is  the  rod  offset  = 19.75  in.  s 1.65  ft. 


Assuming  small  yaw  angles 

Tgt  Fx  i!  (46) 

= 9.25  x 10“7  x 1.65  Y 
= 1.5263  x 10"6  Y 

This  value  was  also  computed  in  the  program  by  computing  yaw  solar  torques 
as  a function  of  yaw  angle.  Using  this  approach  we  obtain 

T , = 1.055  x 10"6  v (47) 

st 

The  value  selected  for  analysis  was  in  between  these  two  values  or 

T = 1.3  x 10“6  y (48) 

st 

The  G(s)  matrix  which  includes  the  effect  of  yaw  solar  torque  is  given  in 
Figure  44. 

Root  loci  were  then  computed  as  a function  of  yaw  solar  torque  magnitude 
for  different  levels  of  control  system  gains.  These  results  are  shown  in 
Figure  45a  through  45d. 
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Figure  44.  G(s)  with  Yaw  Solar  Torque  Effect 

As  seen  in  Figure  45a  for  nominal  yaw  and  roll  control  system  gains,  a 
solar  torque  10  times  larger  than  nominal  would  be  needed  to  cause  the 
oscillatory  pair  to  move  into  the  right  half  plane.  If  the  yaw  gain  was  halved, 
the  solar  torque  would  have  to  be  2 to  4 times  larger  than  nominal.  This  is 
shown  in  Figure  45b.  Halving  the  roll  gain  with  the  yaw  gain  halved  increases 
the  margin  to  10  times  nominal  as  shown  in  Figure  45c.  Finally,  doubling 
the  roll  gain  with  nominal  yaw  gain  would  require  a solar  torque  4 times 
nominal  to  produce  an  instability. 

In  summary,  the  yaw  torque  due  to  solar  pressure  on  the  rods  is  insufficient 
to  generate  an  instability.  The  magnitude  of  the  variables  in  the  force  equa- 
tion cannot  be  varied  realistically  to  produce  no  more  than  a 20  percent  in- 
crease in  solar  force  which  is  far  below  the  200  percent  increase  necessary 
for  instability. 

A solar  force  generated  rolling  moment  due  to  yaw  is  also  present  due  to 
the  imbalance  between  the  upper  and  lower  rods;  however,  it  is  two  orders 
of  magnitude  less  than  the  yaw  moment  and  consequently  insignificant. 
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Figure  45.  Root  Migration  as  a Function  of  Yaw  Solar  Torque  (concluded) 


Magnetic  Torques --Magnetic  torques  are  generated  on  the  spacecraft 
through  the  following  relationship 

T„  = MxB  (50) 

M 

where  M is  the  vector  dipole  moment  of  the  spacecraft  and  B in  the  flux 
density  of  the  earth's  magnetic  field  (centered  tilted  dipole  model).  The 
magnetic  torque  model  has  been  described  previously  and  details  will  be 
omitted  here.  The  flux  density  in  an  earth  reference  frame  is  given  by 

Be  = - (-j^)3  (g 1 1 sin  X - h11  sin  0 ) (51) 

x 

Re  3 

be 

y 

R e 3 o 1 1 

BE  = -2(^")  sin  6 +g!  cos  0 cos  X + hx  cos  0 sin  X) 

z 


(-^■)  (g1°  COS  e + g1  sin  e COS  X + l^1  sin  0 sin  x) 


Choosing  Latitude  ( 0)  = Longitude  (x)  = 0.0  we  obtain 


E 


= 0 


(52) 


R 


BEy  ■ <Ee>3 

be  = '2  <TT>3  <*/> 

z 


Substituting  numerical  values 

B_  = 0 (53) 

E 

x 


R n 

= niQ7Q9/_9  HAH  1 9 v 1H  \ = _A  1 O 9Q  v in  ' 


Now  transforming  from  earth  reference  frame  to  the  local  vertical  and  from 
the  local  vertical  to  the  body  frame  or 


B„  = T 


LV/B  E/LV  E 


sin  63°  - cos  63°  0 

cos  63°  sin  63°  0 


+ 1.9035  x 10 


-4. 1929  x 10 
5.9686  x 10' 


3.7359  x 10 
5 . 9686  x 10  _ 

cos  0 cos  Y cos  0 sin  y 

cos  Y sin  0 sin  0 cos  y cos  0 + sin  y 
- sin  y cos  0 sin  0 sin  0 

cos  y sin  0 cos  0 sin  y sin  0 cos  0 
+ sin  y sin  0 - cos  y sin  0 


-sin  0 
cos  0 sin  0 

cos  0 cos  0 


Assuming  small  angles  and  neglecting  second  order  effects 
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Now  calculating  magnetic  torques  in  the  body  axis 


i 

3 

k 

(58) 

K 

m 

X 

K 

m 

y 

K 

m 

z 

bb 

X 

bb 

y 

bb 

z 

= K (1.9035  x 10_7e  - 3.7359  x 1O_70  + 5.9686  x 10-8) 

B m 

x y 

-K  (-1.9035  x 10"7y  + 3.7359  x 10"7  + 5.9686  x 1O“80) 
m v 

z 

T„  = -K  (1.9035  x 1O"70  - 3.7359  x 1O?0  + 5.9686  x 10"8) 

B m v 

y x 

+K  (1.9035  x 10“7  + 3.7359  x 10"7y  - 5.9686  x 1O~80) 
z 

T = K (-1.9035  x 10"?f  i- 3.7359  x 10“7  + 5.9686  x 1O"80) 

B m v 

z y 

Km  (1.9035  x 10"7  + 3.7359  x 10"?Y  - 5.9C86  x 1O"80) 

my 

where  equals  a scaling  factor  on  the  magnitude  of  the  spacecraft  dipole 
moment. 

Ignoring  the  pitch  component,  G(s)  which  includes  the  effect  of  magnetic 
torques  is  given  in  Figure  46. 

A root  locus  was  computed  as  a function  of  spacecraft  dipole  moment  assum- 
ing Km  = Km  . The  results  are  shown  in  Figure  47.  As  can  be  seen,  a 
spacecraft  dipole  moment  40  times  the  value  of  spacecraft  limitation  is 
required  for  an  instability. 
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Figure  46.  G(s)  Including  the  Effect  of  Magnetic  Torques 


In  summary,  no  stability  problems  are  expected  from  the  effect  of  magnetic 
torques. 


Wheel  Misalignments --As  shown  in  Figure  48,  wheel  alignment  is 


defined  by  a cant  angle  a , and  a wheel  yaw  angle,  Y^y.  which  defines  the 


wheel  position  in  the  satellite  x-y  plane.  The  momentum  of  the  wheels 
in  body  axes  is  given  by 


HW  = _HW  cos  % COSai  ~ HW  cos  % COSa2 
x 1 1 2 2 


(59) 


-HW  cos  % cosa3  - Hw  cos  % cos  “4 
3 3 4 4 


Hyy  ”^W  COS  a 1 ""  COS  a2  ~ COS  ff3  " COS  ff4 

z 1 2 3 4 


The  alignment  specification  requires  alignment  errors  of  no  greater  than 
, o 


1 with  respect  to  a and  Y^..  In  analyzing  the  effect  of  wheel  misalignment, 
it  can  be  shown  that  the  control  equations  reduce  to 


TI  = Kx  (ai^  + a2^  + Kz  + b2Y) 


(60) 


TI  = Kz  (b!^  + b2^  _ Kx  (ai0  + a2Y) 


(61) 


where  the  coefficient  a^,  a^,  b^,  b^,  represented  the  effect  of  wheel 
misalignment.  As  can  be  seen,  these  equations  are  very  similar  to  those 


111 


obtained  analyzing  yaw  sun  sensor  coupling.  The  effect  on  the  a^  coefficient 
of  the  characteristic  equation  is  also  identical. 

Since  the  effect  of  coupling  is  almost  completely  cancelled  out,  the  combina- 
tion of  wheel  misalignment  angles  were  chosen  which  reduced  the  yaw  loop 
gain  and  increased  the  roll  loop  gain.  A root  locus  was  then  computed  as 
a function  of  magnitude  of  misalignment  angles  at  this  configuration.  The 
results  are  shown  in  Figure  49.  As  can  be  seen,  even  a 5°  misalignment 
has  an  insignificant  effect  on  the  location  of  the  oscillatory  pair. 
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3 
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45° 
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4 
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Figure  48.  Wheel  Orientation 
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Figure  49.  Root  Migration  as  a Function  of  Wheel  Misalignment 


IV-B-2.  On  Orbit  Simulation  Results 


Simulation  of  on  orbit  satellite  performance  was  conducted  at  4 sun  line  to 
orbit  plane  angles  (0°,  12°,  45°,  86°)  under  nominal  and  off-nominal  condi- 
tions. Also  simulated  were  independent  failures  of  wheels  1 and  2,  in  both 
compensated  and  uncompensated  modes.  Table  15  gives  definition  of  the 
various  plotted  variables  in  Figures  50  through  75.  Table  16  shows  a log 
of  the  simulation  runs  and  the  input  parameters  for  these  figures. 


Table  15 . Plotting  Symbol  Definitions 


Roll  Angle 
Pitch  Angle 
Yaw  Angle 

Roll  Rate 
Pitch  Rate 
Yaw  Rate 


Roll  Momentum 
Pitch  Momentum 
Yaw  Momentum 


local  vertical  Euler  angles, 
yaw,  pitch,  roll  sequence 


body  rates 


total  satellite  momentum 
in  inertial  space 


Wheel  1 Momentum* 
Wheel  2 Momentum 
Wheel  3 Momentum 
Wheel  4 Momentum 


wheel  momentum  in  body 
space 


Roll  Solar  Torque 
Pitch  Solar  Torque 
Yaw  Solar  Torque 

Panel  Angle 

Sensed  Yaw  Angle 


^ solar  torques  in  body  axis 

( actual  panel  angle 
( output  of  yaw  sun  sensor 


^Replaced  in  later  runs  with  wheel  momentum  by  axis 
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Table  16.  Run  Log  (continued) 


Runs  1,  2,  3,  and  6A  present  satellite  performance  under  nominal  conditions 
at  the  45°,  86°,  0°  and  12°  sun  line  to  orbit  plane  angles  respectively.  As 
shown  in  Figures  50  through  52(  satellite  behavior  is  stable  for  the  first 
three.  There  was  no  calcomp  plot  made  for  6A.  Tracking  errors  are  given 
in  Table  17.  Run  6A  was  made  as  a check  against  the  analytical  results 
which,  while  not  indicating  a problem,  show  reduced  margin  to  yaw  solar 
torques  at  reduced  yaw  gains.  The  12°  sun  line  to  orbit  plane  angle  was 
chosen  for  simulation  since  we  incorrectly  believed  that  the  yaw  gain  in- 
crease switch  was  only  enabled  when  the  sun  line  to  orbit  plane  angle  was 
less  than  10°.  Hence  the  yaw  gain  was  kept  at  its  nominal  value  of  .4  ft -lb/ 
rad.  This  admittedly  is  an  unrealistic  run;  however,  it  does  demonstrate 
that  even  with  very  low  yaw  gains,  no  instability  was  observed.  The  dead 
bands  were  also  removed  since  they  appear  to  add  stability.  As  can  be  seen 
tracking  errors  in  pitch  and  roll  are  comparable  to  the  other  cases.  Yaw 
error  is  larger  which  is  expected. 

The  remaining  runs  were  made  with  off  nominal  conditions  which  will  now 
be  described. 

Earth  Sensor  Misalignments --Two  runs  (1A,  IB)  were  made  with  the  modi- 
fied version  of  the  earth  sensor  as  described  in  the  modeling  section.  Run 
1A  shown  in  Figure  53  is  a r erun  of  Run  1 with  the  only  difference  being  use 
of  approximated  earth  sensor  feedback  in  the  control  system.  Run  IB, 

Figure  54,  contains  a 1°  misalignment  in  the  pitch  and  roll  axes  of  the  earth 
sensor.  Tracking  errors  are  summarized  in  Table  18.  Run  2A  could  also 
be  interpreted  as  perfect  earth  sensor  alignment  but  with  a 1°  misalignment 
of  the  gravity  gradient  rods  in  pitch  and  roll. 
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Table  17.  Tracking  Errors--Nominal  Conditions 


Run  No. 

Sun  Line 

Roll  Error 

Pitch  Error 

Yaw  Sun 
Tracking 
Error 

1 

45° 

.35° 

.37° 

7.0° 

2 

O 

CO 

CO 

. 22° 

.24° 

o 

o 

3 

0° 

. 23° 

.37° 

- 

6A 

12° 

.2° 

.34° 

o 

CO 

o 

Table  18.  Tracking  Errors --Off  Nominal  Conditions 


Run 

No. 

Sun 

Line 

Peak  Roll 
Error 

Peak  Pitch 
Error 

Peak  Yaw 

Sun  Tracking 
Error 

Comments 

1A 

45° 

.40° 

.46° 

7.0° 

with  CES 

IB 

1.66° 

1.04° 

7.0° 

1°  CES  mis- 

4 

45° 

.47° 

.45° 

14.0° 

alignment 

4A 

45° 

.48° 

.425° 

4.47° 

4B 

45° 

.47° 

.47° 

4.8° 

4C 

45° 

.495° 

.476° 

4.8° 

5 

86° 

.306° 

.39° 

3.2° 

6 

0° 

- 

6B 

12° 

.366° 

.585° 

25.0° 

no  deadbands 

6C 

12° 

.334° 

.494° 

18.6° 

no  deadbands 

7 

0° 

Lost  ea 

rth  after  six  c 

rbits 

0° 

.813° 

2.72° 

- 

45° 

.433° 

.804° 

7.5° 

45° 

.52° 

.46° 

3.3° 

86° 

.302° 

.375° 

2.87° 

119 


Yaw  Sun  Sensor  Errors --As  discussed  in  the  modeling  section,  yaw  sun 
sensor  accuracy  is  bounded  by  an  accuracy  envelope  at  25 °C  and  also  by  a 
gain  variation  as  a function  of  temperature.  In  the  simulation  program, 
the  nominal  or  either  boundary  of  the  sensor  curves  may  be  used.  In 
Table  14,  these  curves  have  been  designated  as  curve  1,  2,  or  3 as  follows: 

Curve  1 - lower  than  nominal 

Curve  2 - nominal 

Curve  3 - higher  than  nominal 

In  addition,  the  temperature  variation  can  cause  as  much  as  a ±50  percent 
variation  in  sensor  output.  The  nominal  scale  factor  in  transforming  the 
sensor  output  (microamps)  into  degrees  is  1/6.33.  Hence  1/4.22  repre- 
sents a 50  percent  increase  in  the  sensor  output  and  1/12.66  represents  a 
50  percent  decrease.  Run  4,  Figure  55,  shows  the  effect  of  a 50  percent 
reduction  in  sensor  output.  The  plot  labeled  "sensed  yaw"  is  approximately 
half  the  actual  yaw  tracking  error.  Tracking  error  is  shown  in  Table  18. 
Run  4A,  Figure  56,  shows  the  effect  of  increasing  the  yaw  gain  by  50  per- 
cent. As  shown  in  Table  18,  pitch  and  roll  tracking  error  are  close  to  the 
nominal  performance  for  the  two  runs.  Yaw  performance  varies  propor- 
tionally with  the  gain. 

Run  6B  is  with  the  worst  case  conditions  applied  to  Run  6A  under  the  same 
incorrect  yaw  gain  procedure.  Yaw  gain  is  reduced  by  50  percent  of  its 
nominal  value  over  the  whole  run.  Again  the  dead  bonds  have  been  removed. 
Performance  is  degraded  as  shown  in  Table  18,  but  still  stable.  Run  6C, 
Figure  57,  shows  the  effect  of  the  correct  yaw  gain  implementation.  Per- 
formance is  improved  slightly. 
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Wheel  Misalignments — The  NRL  specifications  permit  a 1°  misalignment  of 
the  reaction  wheels  in  yaw  and  cant  angle.  The  alignment  angles  that  were 
used  in  the  simulation,  as  given  in  Table  15,  were  determined  to  be  worst 
case  combinations  for  the  NDS  satellite.  The  analysis  shown  in  the  previous 
section  indicates  they  may  not  be  the  worst  case  for  NTS-2;  however,  they 
do  effect  a degree  of  nonlinear  coupling  not  accounted  for  in  the  linear 
analysis. 

Run  4B,  Figure  58,  shows  the  results  of  appending  these  misalignments 
onto  the  conditions  of  4A.  As  shown  in  Table  18,  performance  is  not  signif- 
icantly affected. 

Solar  Array  Center  of  Pressure  Shift--A  0.5"  location  uncertainty  was 
assumed  to  exist  in  the  center  of  pressure  on  the  solar  arrays.  The  com- 
bination of  signs  of  the  cp  shift  were  chosen  to  create  a yaw  solar  torque 
and  a further  imbalance  in  pitch  solar  torque.  The  results  are  shown  in 
Run  4C,  Figure  59,  and  Table  18.  Again  performance  is  not  significantly 
affected. 

The  effects  of  a satellite  center  of  mass  variation  ( Ax  = Ay  = . l",  A z = -.  2"), 
a doubling  of  the  spacecraft  magnetic  dipole  moment,  and  a 1°  solar  array 
tracking  error  were  also  included  in  Runs  4,  4A,  4B,  and  4C.  The  Az  sign 
was  chosen  to  increase  the  pitch  solar  torque  imbalance.  The  dipole  moment 
variation  is  implemented  through  a scale  factor  parameter  for  changing 
Newton-meters  to  ft -lbs.  For  a dipole  moment  of  1,000  pole  cm,  the  scale 
factor  is  .73.  Since  the  one  body  simulation  used  in  orbit  runs  assumes 
perfect  panel -tracking,  a 1°  tracking  error  was  inserted  to  simulate  the  1° 
dead  band  in  the  solar  array  drive.  All  of  these  perturbations  contributed 
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to  the  slight  degradation  in  tracking  performance  as  shown  in  Table  18.  The 
remaining  runs  contained  all  the  above  described  perturbations  for  generating 
"worst  case"  conditions.  Run  5,  Figure  60,  is  a repeat  of  Run  4C  at  the  86° 
sun  line  condition.  Slightly  better  performance  is  shown  in  Table  18.  Run 
6,  Figure  61,  shifts  the  simulation  to  noon  turn  and  puts  the  yaw  gain  to  50 
percent  of  its  nominal  value. 

The  final  five  runs  simulate  the  result  of  wheel  failures,  both  compensated 
and  uncompensated  under  worst  case  condition.  Run  7 has  an  uncompensated 
wheel  1 failure  at  the  noon  turn  condition.  Only  line  printer  plots  are  avail- 
able for  this  run  and  have  not  been  included.  The  satellite  did  lose  the  earth 
at  approximately  6 orbit  into  the  run.  Run  7A,  Figure  62,  presents  the  same 
condition  as  Run  7 however,  with  the  compensation  logic  included.  Perfor- 
mance does  degrade  seriously  as  shown  in  Table  18,  The  large  attitude 
excursion  in  pitch  at  one  point  in  the  run  is  felt  to  be  due  to  an  instantaneous 
wheel  unloading.  Wheel  2 failures  at  two  less  demanding  mission  conditions 
were  also  simulated.  Run  8,  Figure  63,  shows  uncompensated  wheel  2 
failure  response.  Performance  is  surprisedly  good.  Run  9,  Figure  64, 
shows  the  effect  of  adding  compensation.  Performance  does  improve, 
particularly  in  yaw  tracking.  Run  10,  Figure  65,  is  identical  to  Run  9 
except  at  the  86°  sun  line  condition.  Again  performance  is  slightly  better. 

Summary  --There  appears  to  be  no  problem  in  orbit  performance  of  the 
satellite.  Tracking  error  is  generally  less  than  .5°  in  pitch  and  roll  and 
stable  in  yaw.  An  uncompensated  wheel  1 failure  will  lead  to  the  earth 
falling  out  of  the  field  of  view  of  the  earth  sensor  if  not  compensated  for  in 
6 orbits.  Compensated  wheel  failure  performance  under  worst  case  condi- 
tions indicates  transient  attitude  errors  of  2.70°.  These  do  not  appear  to 
be  a problem. 


FV-C.  DEPLOYMENT  AND  ACQUISITION 


In  determining  permissible  initial  attitudes  and  rates  from  which  acquisition 
can  be  achieved  after  deployment  of  the  rods,  it  is  necessary  to  specify  an 
acquisition  requirement.  One  reasonable  requirement  is  that  the  initial 
rates  and  attitudes  be  such  that  acquisition  be  achieved  without  saturating 
the  earth  sensor  and  yaw  sun  sensor.  If  this  requirement  is  to  be  met,  the 
satellite  wheels  must  be  able  to  absorb  the  momentum  in  the  satellite  due  to 
attitude  rates.  In  addition,  the  wheels  must  provide  the  satellites  with  the 
orbital  rate  components  necessary  for  tracking  the  earth,  and  the  sun  track- 
ing rates  for  tracking  the  sun.  This  latter  requirement  is  a function  of  sun 
line  position  and  consequently  the  yaw  angle  necessary  for  sun  pointing.  If 
the  yaw  angle  is  0 , the  orbital  rate  component  must  be  built  up  in  the  satel- 
lite pitch  axis.  If  the  yaw  angle  is  90°,  then  the  satellite  roll  axis  will  see 
the  effect  of  orbital  rate.  Between  0°  and  90°  yaw  angles,  it  is  a combina- 
tion of  the  two.  The  momentum  transfer  procedure  may  be  analyzed  as 
follows.  The  ACE  system  commands  wheel  momentum  as  a function  of 
attitude.  Hence  the  maximum  wheel  momentum  that  can  be  commanded 
will  be  when  the  sensors  are  at  the  saturation  point,  or 

Hwx  = Kx  x 4°  = 3,75  x 5^r~3  = ft -lb -sec 
max  ' * 


' ■ K *4°  =3.75 

ymax  Y 

rz  x30°-  0-4*rr3 

max 


, 26  ft  -lb-sec 


= . 21  ft -lb-sec 


These  relationships  assume  no  initial  attitude  error. 


Allowable  rates  for  achieving  acquisition  can  then  be  computed  from 
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If  the  acquisition  procedure  was  not  to  be  initiated  until  after  the  rods  were 
fully  deployed,  the  allowable  for  nominal  conditions  would  then  be 

Vax  ^°St  deployment)  = .0075°/sec 

Vax  ^°St  deployment)  = .0078°/sec 

rmax  ^Post  deployment)  = . 06°/sec 

Since  orbit  rate  is  . 008°/sec,  these  results  indicate  that  it  would  be  diffi- 
cult at  yaw  angles  of  0°  or  90°  to  acquire  without  saturating  the  earth  sensor 
and  that  assuming  no  initial  body  rates. 

In  the  runs  we  simulated,  we  engaged  wheel  control  at  the  start  of  deploying 
the  rods.  At  that  condition  the  allowable  rates  are 

Pmax  ^Pre  deployment)  = ,09° /sec 

qmax  (pre  deployment)  = .191° /sec 

r (pre  deployment)  = ,06°/sec 

Hence,  if  the  wheels  are  to  be  turned  on  at  the  same  time  as  the  rods  begin 
deployment,  the  permissible  rates  will  fall  somewhere  in  between  the  two 
sets  of  values.  (Note:  permissible  yaw  rate  is  the  same.)  This  result  was 
verified  in  simulation  and  will  now  be  discussed. 
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Five  runs  were  made  to  investigate  rod  deployment  and  acquisition  condi- 
tions. These  are  listed  with  their  pertinent  input  parameters  in  Table  15. 
The  86°  sun  line  was  chosen  since  the  orbital  rate  component  appears  almost 
entirely  in  the  roll  axis  at  this  condition.  Run  19,  Figure  66,  and  Run  20, 
Figure  67,  show  acquisition  problems  for  initial  yaw  rates  of  . l°/sec  and 
,05°/sec  respectively.  The  ,05°/sec  rate  is  less  than  the  maximum  per- 
missible. However,  with  a 10°  initial  offset  in  yaw  only  20°  of  the  30° 
sensor  range  was  available.  For  that  condition,  the  permissible  yaw  rate 
is  on  the  order  of  ,04°/sec.  Run  21,  Figure  68,  did  acquire  without  sensor 
saturation.  The  roll  attitude  response  shows  the  effect  of  the  additional 
orbit  rate  term.  Roll  attitude  approaches  4°  whereas  pitch  attitude  reaches 
approximately  1.7°.  When  orbit  rate  is  added  to  the  initial  roll  body  rate, 
the  equivalent  initial  roll  rate  would  be  .015° /sec. 

Run  23,  Figure  69,  shows  the  effect  of  varying  the  rod  deployment  rate  +10 
percent  for  rod  1 and  -10  percent  for  rod  2.  There  is  no  measurable  effect. 

i 

Run  24,  Figure  70,  adds  a 50  percent  reduction  in  yaw  gain  to  the  conditions 
of  Run  23.  Acquisition  is  still  achieved. 

In  summary,  it  has  been  demonstrated  that  initial  conditions  of 


p = .007°/ sec 

body  rates 

q = .007°/sec 
r = .01°/ sec 

Euler 

<t>  = .25° 

angles 

0 = .25° 

Sun  mis- 
alignment 

o 

= .10 
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will  permit  deployment  and  acquisition  withou!  senior  saturation.  Margins 
for  acquisition  at  different  sun  line  conditions  have  also  been  provided. 


IV-D.  ROD  AND  PANEL  FLEXIBILITY  EFFECTS 

The  effects  of  rod  and  panel  flexibility  were  investigated  by  1)  determining 
three  axis  gain  and  phase  margins  via  frequency  response  analysis,  and 
2)  computer  simulation  which  included  the  actual  ACE  system  and  the  com- 
plete three  body  dynamics  model. 

Gain  and  phase  plots  for  the  roll,  pitch  and  yaw  axes  are  given  in  Figures 
76,  77,  and  78  respectively.  The  first  and  second  unsymmetric  rod  bending 
modes  are  visible  in  both  the  roll  and  pitch  plots,  however,  because  of  the 
fifth  order  filtering,  their  magnitude  is  very  small.  (Note-  a structural 
damping  ratio  of  .0005  sec  * was  assumed  in  this  analysis.)  No  elastic 
modes  were  observed  in  the  yaw  axis  because  of  the  small  magnitude  of  the 
yaw  excitation  moment.  A summary  of  the  gain  and  phase  margin  charac- 
teristics is  given  in  Table  19. 

r ive  computer  runs  were  made  which  by  varying  the  run  length  investigated 
both  the  effect  of  torsional  flexibility  of  the  solar  panels  and  also  the  bending 
effects  of  the  rods.  These  runs  are  described  in  Table  20  and  the  response 
histories  are  given  in  Figures  71  through  75.  Symbol  definition  is  given  in 
Table  21  . 

Run  25,  Figure  71.  simulated  two  seconds  of  satellite  performance  with  the 
most  complete  dynamic  model.  No  instabilities  were  observed.  Because 
the  run  was  so  short,  it  is  possible  to  see  the  8 Hz  panel  t-rsional  mode 
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dynamics  both  in  the  modal  c oordinate  plot  itself  and  also  in  the  panel  motion 
at  the  drive  motor  and  at  the  pitch  sun  sensor.  This  run  was  made  without 
the  GE  modification  in  the  solar  array  drive  and  the  motor  torque  pulses  can 
be  seen  at  their  maximum  2 pulse  per  second  frequency.  The  torsional  mode 
is  excited  at  each  pulse;  however,  damping  can  be  observed  between  pulses. 

Runs  26  and  27,  Figures  72  and  73,  are  extensions  in  time  of  Run  25  to  60 
seconds.  Run  27  includes  the  GE  modification  in  the  solar  array  drive  loop. 
No  instabilities  were  observed  in  either  simulation.  It  is  interesting  to  note 
that  even  after  6°  seconds  running  the  wheels  have  still  not  come  on.  This 
is  because  of  the  deadbands  in  the  wheel  loop  and  the  fifth  order  filters  which 
have  equivalent  first  order  time  interval  of  six  minutes.  Attitude  response 
shows  the  effect  of  orbital  rate  which  because  of  the  86°  sun  line  condition 
appears  almost  entirely  in  roll.  Run  26  contains  a plot  of  sensed  angle  error 
showing  that  the  10°  error  has  been  reduced  within  the  1°  deadband  in  approx- 
imately 40  seconds.  Run  27  contains  a similar  plot  except  that  the  output  of 
the  20  second  filter  is  plotted.  Response  is  much  smoother  as  can  be  seen. 
Because  of  the  lag  and  the  pulse  frequency  reduction,  panel  error  was  still 
1.96°  at  60  seconds. 

Run  27A,  Figure  74,  is  a repeat  of  run  26,  assuming,  however,  that  the 
solar  array  is  grounded  at  the  drive  motor  instead  of  free  wheeling.  This 
implies  that  the  rigid  body  can  now  be  excited  by  the  panel  torsion  motion 
which  was  not  possible  in  the  free  wheeling  model.  As  discussed  in  the 
modeling  section,  the  frequency  of  the  first  symmetric  torsion  mode  is 
reduced  from  8 Hz  to  2.6  Hz  for  this  model.  There  is  very  little  difference 
between  the  two  runs.  Rigid  body  pitch  attitude  is  just  slightly  larger.  In 
runs  26,  27  and  27A,  only  one  period  of  rod  bending  effects  are  observable. 
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Tip  mass  deflections  in  the  x and  y directions  are  plotted  and  show  peak 
magnitude  on  the  order  of  . 1 - . 15  feet.  In  order  to  see  the  effects  of  the 
rod  bending,  Run  28,  Figure  75,  was  made.  In  this  run  the  solar  array 
torsion  mode  was  deleted;  however,  the  solar  array  drive  model  was  still 
retained.  The  GE  modification  was  included  in  this  run. 

The  effect  of  the  first  antisymmetric  rod  bending  mode  at  . 12  rad /sec  can 
be  observed  on  the  rigid  body  pitch  attitude  response.  Again  no  instabilities 
were  observed.  Wheel  activity  was  apparent  at  about  70  seconds  into  the 
run.  One  interesting  point  is  that  roll  attitude  has  reached  3.74°  in  600 
seconds  and  appears  to  be  growing.  This  is  in  direct  agreement  with  the 
results  previously  presented  on  deployment  which  stated  that  at  post  deploy- 
ment inertias  the  maximum  permissible  rate  is  ,007°/sec.  Orbital  rate 
which  in  this  case  appears  almost  entirely  in  the  roll  axis  is  ,008°/sec. 
Observation  of  the  roll  rate  plot  will  show  that  roll  rate  has  very  nearly 
reached  the  ,008°/sec  required  which  is  consistent  with  the  500  second 
time  constant  of  the  roll  rate  dynamics.  Hence  roll  attitude  should  be 
flattening  and  gradually  through  gravity  gradient  torques  return  to  zero. 
These  results  indicate  that  despite  using  a modified  version  of  the  control 
system  in  the  deployment  simulation  analysis,  the  results  are  reasonable. 

In  summary,  no  problems  were  detected  in  the  analysis  of  rod  and  panel 
flexibility  effects.  There  are  very  large  stability  (gain  and  phase)  margins 
in  the  design  and  the  simulation  could  detect  no  instabilities  due  to  nonlinear 
eftects. 


Table  19.  ACE  System  Stability  Margins 
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Table  20.  Flexibility  Effects  Run  Table 
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Table  21.  Hexibility  Run  Symbol  Definition 


Roll  angle 

Pitch  angle 

Yaw  angle 

Euler  angles  in  the 
local  vertical 

Roll  rate 

Pitch  rate 

Yaw  rate 

body  rates 

Sensed  yaw  angle 

output  of  yaw  sun  sensor 

First  symmetric 
panel 

first  symmetric  torsional 
mode  in  modal  coordinates 

Theta  at  DM 

position  of  the  solar  array 
shaft  at  the  drive  motor 

UAX  tip  mass 

deflection  of  tip  mass  in 
the  x direction 

UAY  tip  mass 

deflection  of  tip  mass  in 
the  y direction 

Theta  at  sensor 

position  of  the  solar  array 
shaft  at  the  pitch  sun  sensor 

Panel  angle  error 

pitch  sun  sensor  output 

Panel  angle  error 
sensed  lag 

pitch  sun  sensor  output 
filtered  through  Ge  lag 

Figure  52.  Run  3 - Sunline  0°  - Nominal  (continued) 
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Figure  53.  Run  1A  - Sensed  Pitch  and  Roll  Angles 
in  Feed  Back  Nominal 
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Figure  53.  Run  1A  - Sensed  Pitch  and  Roll  Angles 
in  Feed  Back  Nominal  (continued) 
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Figure  53.  Run  1A  - Sensed  Pitch  and  Roll  Angles 
in  Feed  Back  Nominal  (concluded) 
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Figure  54.  Run  IB  - Sensed  Pitch  and  Roll  Angles 
in  Feed  Back  Earth  Sensor  Misaligned 
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Figure  54.  Run  IB  - Sensed  Pitch  and  Roll  Angles  in  Feed 
Back  Earth  Sensor  Misaligned  (continued) 
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Figure  54.  Run  IB  - Sensed  Pitch  and  Roll  Angles  in  Feed 
Back  Earth  Sensor  Misaligned  (continued) 
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Figure  54.  Run  IB  - Sensed  Pitch  and  Roll  Angles  in  Feed 
Back  Earth  Sensor  Misaligned  (continued) 
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Figure  58.  Run  4B  - Sunline  45  - 4A  + Wheel 

Misalignments  (continued) 
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Figure  58.  Run  4B  - Sunline  45  - 4A  + Wheel 

Misalignments  (continued) 


Figure  58.  Run  4B  - Sunline  45°  - 4A  + Wheel 
Misalignments  (continued) 


Figure  59.  Run  4C  - Sunline  45~  - 4B  + CP  Shift  (continued) 


"W 


3 CM  M1E  DCC/XCJ 


198 


Figure  62.  Run  7A  Wheel  1 Failure  Compensated  - Worst 
Case  - ‘ online  0°  (continued) 
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Figure  63,  Run  8 - Wheel  2 Failure  Uncompensated  - Worst 


Case  - Sunline  45° 
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Figure  64.  Run  9 - Wheel  2 Failure  Compensated 
Worst  Case  - Sunline  45°  (continued) 
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Figure  64.  Run  9 - Wheel  2 Failure  Compensated 
Worst  Case  - Sunline  45°  (continued) 
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Figure  64.  Run  9 - Wheel  2 Failure  Compensated  - 
Worst  Case  - Sunline  45  (concluded) 
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Figure  65.  Run  10  - Wheel  2 Failure  Compensated  - 
Wor^t  Case  - Sunline  86°  (continued) 
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Figure  65.  Run  10  - Wheel  2 Failure  Compensated  - 
Worst  Case  - Sunline  86°  (continued) 
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Figure  65.  Run  10  - Wheel  2 Failure  Compensated  - 
Worst  Case  - Sunline  86  (concluded) 
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Figure  68.  Run  2]  - Deploy  (P.  Q.  R.  )-(.  007.  . 007. 01) 
(PHI.  TUT.  P91)  = (.25.  .25. 10.0) 


Figure  68.  Run  21  - Deploy  (P.  Q.  R ) = ( . 007 . . 007 . . 01) 
(PHI.THT.  P91)  = (.25.  .25.  10.0)  (continued) 

243 


Figure  68.  Run  21  - Deploy  (P.  Q.  R.  ) = (.007.  .007.  .01) 
(PHI.THT.  P91)  = (.25..  25.  10.0)  (continued) 


Figure  69.  Hun  23  - Deplo  (P.Q.  R.  ) 
(.25. .25.  10.0)  1,'ODl  and 
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Figure  69.  Run  23  - Deploy  (P.  Q.  R.  ) 
(.  25.  . 25.  10. 0)  ROD1  and 
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Figure  69.  Run  23  - Deploy  (P.  Q.  R.  ) = ( . 007.  . 007 . . 01)  (PHI.THT.  PSI)  = 
( . 25 . . 25 . 10. 0)  ROD1  and  ROD2  Rates  Unequal  (continued) 
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Figure  70.  Run  24  - Deploy  (P.Q.R.  ) = (.007.  .007.  .01)  (PHI.THT.  PSI)  - 


(.25.  .25.10.0)  ROD1  and  ROD2  Rates  Unequal  KSUB(PSI)  = 
KSUB(PSI)/ 2 (continued) 
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Figure  70.  Run  24  - Deploy  (P.  Q.R.  ) = (.007.  .007.  .01)  (PHI.  THT.  PSI) 
(.  25.  . 25.  10.0)  ROD1  and  ROD2  Rates  Unequal  KSUB(PSI)  = 
KSUB( PSI) / 2 (continued) 
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Figure  70.  Run  24  - Deploy  (P.  Q.  R.  ) = (.007.  . 007.  .01)  (PHI.  THT.  PSI)  = 
(.25.  .25.  10.0)  ROD1  and  ROD2  Rates  Unequal  KSUB(PSI)  = 
KSUB(PSI)/2  (concluded) 
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Figure  71.  Run  25-2  Second  Run  3 Body  (concluded) 


Figure  72.  Run  26  - 60  Second  Run  3 Body 
without  C.E.  Fix  (continued) 
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Figure  72.  Run  26  - 60  Second  Run  3 Body 
without  C.  E.  Fix  (continued) 
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Figure  72,  Run  26  - 60  Second  Run  3 Body 
without  C.  E.  Fix  (concluded) 
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Figure  73.  Run  27  - 60  Second  Run  3 Body 
with  C.  E.  Fix  (continued) 
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Figure  73  Run  27  - 60  Second  Run  3 Pody 
with  C.  E.  Fix  (concluded) 


Figure  74.  Run  27A  - 60  Sec  3 Body  Run  with 
Locked  Up  Drive  Motor  (continued) 


Figure  74.  Run  27A  - 60  Sec  3 Body  Run  with 
Locked  Up  Drive  Motor  (continued) 
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Yaw  Wheel  Frequency  Response 


SECTION  V 
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CONCLUSIONS 


To  reiterate  the  comments  given  in  Section  I,  there  appears  to  be  no  problem 
in  the  performance  of  the  NTS-2  system.  The  only  recommendation  we  have 
is  to  consider  the  effect  of  yaw  sensor  uncertainty  and  orbital  rate  considera- 
tions on  deployment  and  acquisition  procedures.  In  conclusion,  we  would 
like  to  compliment  General  Electric  on  an  excellent  job  of  analysis  and  docu- 
mentation of  results. 
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